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RADHESHYAM KHANDELWAL 
Ph.D. 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 
JANUARY, I960 

SOME PROBLEMS IN THE SLOSHING OF 
LIQUID IN MOVING containers 


The sloshing of liquids in partially filled moving 
containers occurs in a diverse class of problems. The moving 
container may be a fuel tank of an aircraft, a launch vehicle 
or an automobile} a ground based vehicle, such as, a truck, 
a railway wagon or an ocean going ship carrying liquid as 
cargo; an overhead tank, lake or dam, during ground motion, 
due to an earthquake ^ The container motion causes the 
contained liquid to oscillate. The oscillating liquid causes 
time dependent forces and moments to act on the container 
walls. Further, if the container walls are flexible, then 
there is an interaction between liquid oscillation and 
vibration of the container walls.- For the structural design 
of the container, it is necessary to know the time-history 
of the forces and moments acting on the container walls due 
to sloshing. In the case of vehicles carrying liquid, the 
sloshing may be induced due to acceleration or deceleration 
of the vehicle and/or, the vibration of the vehicle. The 
sloshing of liquid generates time dependent forces and moments, 
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which in turn affect the vehicle vibration. Thus sloshing 
is coupled with vehicle vibration and forms an integral psrt 
of the overall vehicle dynamics. Under certain conditions 
the liquid in a container may undergo increasingly large 
amplitude oscillations leading to an interesting class of 
problems of instability. 

For a mathematical treatment of the sloshing of liquid 
in a container, it is necessary to derive the equations govern- 
ing the motion of the container and the liquid. The governing 
equations consist of a system of coupled equations of hydro- 
dynamics and the equations of vibration of the container system. 
The exact nature of equations depends upon the physical problem. 
Generally, however, the equations are coupled, non-linear, 
partial and ordinary differential equations . Under certain 
assumptions, the equations of hydrodynamics can be linearized. 
Even these linearized equations are not always amenable to 
analytical treatment. For a large class of problems, parti- 
cularly those concerning vehicle dynamics, it is convenient 
to replace the oscillating liquid by an equivalent discrete 
element mechanical model. For such a. model, the governing 
equations consist of a system of coupled ordinary differential 
equations. These equations again may be nonlinear and may be 
further linearized or solved numeri cally »| 

The sloshing of liquid in containers mounted on ground 
based vehicles, constitutes a special class of problems. 

In these problems, the basic sources of vibration are the 
acceleration or ueceleration of the vehicle and/or, the 
ground induced excitation due to the unevenness of the 
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prepared surface, such as, road, railway track or runway. 

Due to random nature of the unevenness of the prepared 
surfaces on which the vehicle rolls, the ground induced 
excitation is random in nature, 'i'he unevenness of prepared 
surfaces can be treated as a homogeneous one -dimensional 
random process with the space co-ordinate along the direction 
of motion as the indexing parameter. If the vehicle moves 
with a constant velocity, the base induced excitation is 
stationary and the vehicle response may be determined using 
the well known techniques of random vibration. If the vehicle 
moves with variable velocity, the base induced excitation 
becomes non-st ationary . An analytical treatment of the non- 
stationary problem is much more difficult. Since the under- 
lying process is homogeneous in spac-e\ co-ordinates , the 
concept of evolutionary random process can be applied by 
transforming the governing equations from time to space co- 
ordinates. On such a transformation, the differential equations 
have variable co-eff i cients , but the excitation is homogeneous. 
!he equations can now be solved either analytically in the 
frequency domain, using the concept of evolutionary spectra, 
or through simulation. 

I'he present investigation deals with following problems 
in the sloshing of liquid in moving containers: 

(i) Investigation of the stability of the plane free surface 
of liquid in a. rectangular container with flexible 
bottom, under periodic vertical excitation. 
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(ii) Development of a mechanical model for liquid sloshing 
in a rectangular container under the simultaneous 
action of vertical excitation and horizontal motion. 

(iii) Study of the dynamics of a vehicle carrying liquid as 
cargo and comparison of the results with a vehicle 
carrying an equivalent solidified mass of the liquid 
under identical random track induced excitation. 

An incompressible inviscid liquid in a. rectangular 
container with flexible bottom is subjected to longitudinal 
periodic excitation, liquid motion is assumed to be linear. 
The container bottom is assumed to be simply supported plate 
on all edges and undergoes elastic deformations under thin 
plate approximations. The coupled differential equations have 
been solved to determine the regions of parametric instabi- 
lities. Numerical results have been presented to determine 
the primary and secondary instability regions for different 
cases of natural frequencies of liquid and bottom wall of the 
container and compared with stability regions for a container 
with rigid bottom. 

An equivalent mechanical model is constructed for the 
liquid sloshing in a rigid rectangular container. The 
container is subjected to a constant horizontal acceleration 
and a periodic longitudinal acceleration. The liquid is 
assumed to be incompressible and inviscid. The governing 
hydrodynamic equations are linearized. The equivalent 
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mechanical model is determined by the equivalence of forces 
and moments developed by the sloshing liquid on the walls of 
the container. I'he mechanical model has a rigid mass and 
infinite number of oscillating pendulum masses. I'he analogy 
is valid for stable free surface conditions of the liquid. 

In the last phase of the investigation, the dynamic 
response of a railway wagon carrying liquid as a cargo is 
studied. Two models of the wagon are considered. In the 
first the wagon is assumed to respond under the heave mode 
and in the second it responds under the heave and pitch modes. 
Two types of motions of the wagon are considered; (a) motion 
with constant acceleration, and (b) motion with constant 
velocity. 

The equations governing the vibration of the wagon 
models are nonlinear. These equations are solved for the 
heave model moving with constant acceleration, through 
simulation. The track unevenness is simulated on a digital 
computer assuming a suitable power spectral density of the 
track. The liquid cargo is replaced by an equivalent 
mechanical model as described earlier < The equations govern- 
ing the vehicle-fluid system are coupled and nonlinear. The 
response is determined by numerical integration in time 
domain for different heights of the liquid in the container. 

It is found that the vehicle response is a non-st at ionary 
narrow band random process. The liquid response reaches a 
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steady state quite early and the relative wave height, at 
the -free surface of the liquid in the container, decreases 
with increase in the depth of the liquid- The coupling 
effect is found to he negligible. 

The response of the wagon is also determined for a 
model having heave and pitch degrees of freedom. The liquid 
cargo is replaced by its equivalent mechanical analog. 

Equations of motion of the vehicle-fluid system are nonlinear 
and coupled. They contain inertial, velocity and restoring 
force nonlinearities. Equations of motion are linearized by 
one of the methods of stochastic linearization. The linearized 
equations of motion are transformed in space domain. The 
base excitation is assumed to be homogeneous in space domain. 
Assuming a suitable power spectral density of the track 
unevenness , the solution is obtained in the frequency domain. 
The response statistics are determined and compared with 
the vehicle system which has an equivalent solidified mass of 
the liquid. 

The wagon is assumed to move with constant velocity and 
constant acceleration. For constant -velocity case, the 
results obtained through the analytical form of the response 
statistics are presented. For constant acceleration case the 
equations of motion are solved numerically using the simulation 
technique . 
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When the wagon moves with constant velocity, the vehicle 
response statistics vary rapidly in the neighbourhood of small 
vehicle velocities. The response settles down with increase in 
vehicle velocity. it large vehicle velocities the response of 
-the vehicle, in hea.ve and pitch modes, tends to be independent 
of vehicle velocity. The relative wave height of the liquid -at 
the free surface is found to increase with decrease in depth 
of the contained liquid. The relative wave height of the 
liquid increases with the increase in the eccentricity of the 
container c.g. placed with respect go vehicle c.g. The same 
is true for the shear force and bending moment developed by the 
sloshing liquid at the bottom edge of the side wall of the 
container. 

On comparing the response of vehicle-liquid system with 
that of the vehicle-solidified mass system, it is found that 
the sloshing of liquid has the effect of reducing the vehicle 
response over most of the velocity range of interest. At some 
velocities, sloshing of liquid causes increase in vehicle 
response, but the range of such velocities is narrow and occurs 
at small vehicle velocities. The decrease in vehicle response 
is due to absorber— effect . The liquid absorbs the energy from 
the vehicle thus reducing the response, which in turn affects 
the liquid sloshing. The absorber-effect increases with 
increase in eccentricity of the container c.g. placed with 
respect to vehicle c.g. For a vehicle moving with constant 
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acceleration the overall vehicle and liquid response is simi- 
lar to the constant velocity case. The reduction in vehicle 
response due to transfer of energy to the liquid is significant 
at high vehicle velocities. It. is concluded that neglecting 
liquid sloshing provides conservative estimates of forces and 
displacements for under carriage design. 



CHAPTER 1 


INTRODUCTION 

1 • 1 ‘£M Sloshing of Liquids 

I'he dynamic behavior of liquids in moving containers, 
also known as sloshing, is of considerable interest to 
engineers, scientists, geophysicists, seismologists and 
other scientific workers . The spilling of a partially 
filled container carried hurriedly from one place to another 
is a common experience. Drinking a cup of tea or coffee 
can be a frustrating experience while riding on a vehicle 
moving over an uneven surface. Similar phenomena of liquid 
sloshing are observed in fuel tanks of aircrafts, rockets, 
automotive vehicles, large ocean going ships and rail-road 
tanm cars. The motion of water in wells, lakes and dams 
during earthquakes constitutes large scale sloshing of a 
confined liquid mass. 

The phenomena of liquid motion within a. containing 
vessel presented a challenge to ohe hydro dyn ami cists of the 
past century. Their interest was stimulated by the need 
to study the earthquake effects on large dams, gasoline 
oil tanks, water reservoirs and elevated water tanks 
[31,49,50]. In applying the terrestrial molten core theory, 



10 


geophysicists have been interested in The study of motion 
of liquids in rotating containers [33 3 • the design of 
Boiling— Wat er- React or (BwR) [ 3 ] requires that the water 
level does not drop below a certain level under earthquake 
excitation, ihis is primarily due to The danger of escaping 
superheated steam through the iuret pipe. Similarly; the 
dynamic response analysis of a fast breader reactor should 
consider the sloshing of liquid sodium* used as coolant, in 
containers [9]. 

i'he problem of liquid, fuel response in containers of 
aircrafts, rockets and missiles and its effect on vehicle 
dynamics has received considerable attention during the 
last three decades. Dynamic stability, flutter, control- 
lability, attitude stabilization, elastic deformations 
produced due to forces created by fuel are some of the 
problems of continuing interest in aerospace applications. 

Ground based vehicles often carry liquid as fuel or 
cargo. these include the rockets during their motion on 
a sled, aircrafts during taxi, take-orf and landing-, 
automotive vehicles and. railway wagons. these vehicles 
are subjected mainly to two types of excitations, the 
first is the vertical base excitation of the vehicle 
Through the points of contact on an uneven surface on 
which the vehicle is moving. i'ho second is the 
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horizontal acceleration due to varying velocity of the 
vehicle. The dynamic analysis of such vehicles is consi- 
dered primarily to understand the design aspects of the 
vehicle structure, passenger comfort, cargo safety and 
contro liability of the vehicle. 

As pointed out earlier, the liquid in a container 
mounted on a vehicle is subjected to a vertical excitation 
due to surface unevenness. This causes parametric excita- 
tion of the fluid mass and. under certain conditions the 

Yt 

free surface may undergo violejf oscillations. This may 

5 

result in disintegration of the surface and entrainment of 
small vapor bubbles may cause choking of the pipes connected 
to the container. Thus a study of the stability of the 
free surface of fluid in a container is also of considerable 
interest [1-3,8] . 

1.1.1 The design and dynamics of vehicles 

The modern trend in vehicle design is towards higher 
speed, increased structural flexibility, light weight, low 
cost and longer life. The design of vehicle structure is 
strongly influenced by the dynamic loa.ds due to vehicle 
motion. Structural components of vehicles, such as, landing 
gear and wing root in aircrafts, are prone to fracture or 
fatigue failure due to dynamic loading. 
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A vehicle should offer, psychological and physical 
comfort and safety to the passengers. 1'he induced accelera- 
tion and the range of frequency are the two important 
parameters affecting the safety and comfort. Certain 
cargo, such as fragile materials, cannot withstand vibration 
effects beyond, a certain level. In such a situation, 
operating vehicle speed, becomes an important design cons- 
traint and this in turn influences the design of the under- 
carriage. Controllability and manoeuvrability of the 
vehicle by the driver are other important design aspects. 
Generally, gradual loss of control is followed by increasing 
velocity. Similarly, the manoeuvrability is reduced 
whenever there is a loss of ground contact or contact 
pressure of the wheels . 

Ground based, vehicles usually roll over an uneven 
surface and are subjected to base excitation. !The charac- 
t eristics of base excitations are determined by the type 
of contact, ground unevenness and. vehicle speed, hue to 
random nature of ground unevenness, a venicle moving on 
such a surface is subjected to random excitation. It is, 
therefore, necessary to carry out a random vibration analysis 
to determine the dynamic response of vehicles. 

I he vehicle dynamic response is determined by the 
nature of base excitation and the dynamic characteristics 
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of the vehicle, primarily, its natural frequency, modeshapes 
and damping. In vehicles carrying fluid in containers, 
vehicle motion and ground induced excitations cause fluid, 
to oscillate* which in turn causes time varying loads to 
act on the vehicle, Thus, there is a coupling between 
liquid sloshing and vehicle response. Under certain condi- 
tions, this interaction may be significant both from the 
point of view of vehicle dynamics and liquid sloshing. 

1 . 2 Review of Previous Work 

The present work deals with some aspects of liquid 
sloshing and, also its effects on vehicle dynamics. A 
brief review of previous work in liquid sloshing and vehicle 
dynamics is given in the following sections. 

1.2.1 Hydrodynamic effects induced in containers 

i'he governing equations of motion for the liquid 
sloshing are generally non-linear in nature [1,2,8,17]. 
However, under certain approximations tnese equations can be 
linearized [1,8,31]. There are several interesting aspects 
of liquid, sloshing, such as, liquid behavior under gravita- 
tional force, the bubble formation and. growth, fluid-elastic 
interaction, vortex formation in draining containers, 
spinning containers, connected containers and other related 
problems. Cooper [17] and Abramson [2] have given a detailed 
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review of dynamic behavior of liquids in moving containers. 

A comprehensive review and analysis is given in Ref. [1] . 

In the sequel, a brief review of linearized liquid behavior 
in a container is given under the gravitational force. 

1 . 2 . 1 . 1 Linear inviscid fluid oscillations in rigid containers 

Scientists [31,41,8], seismologists [24], aerospace 
engineers [22,29,36] have shown considerable interest in 
developing comprehensive linearized theory of fluid oscilla- 
tions in containers. For an incompressible fluid, Graham et al . 
[22] gave a detailed analysis based on horizontal, translation, 
pitching about a horizontal axis and rolling about a vertical 
axis for a partially filled rectangular container. Similar 
analysis for cylindrical containers is reported in [1-3,17, 

29 ] . Workers in reservoir and reservoir-dam systems also 

have shown keen interest to take into account the hydrodynamic 
effects on dams during earthquake excitations for incompres- 
sible [49] and compressible [15,38] nature of water. 

1 . 2 . 1 . 2 Equivalent mechanical models 

To study the dynamic forces coming on the liquid-filled 
structures, it is convenient to replace the liquid by an 
equivalent mechanical system. The reason is that the discrete 
mass, stiffness and damping elements are more conveniently 
included in the equations of motion for the entire system. 
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Generally, an equivalent mechanical system is constructed 
by the equivalence of forces and moments developed by the 
fluid on the walls of the container. Such mechanical systems 
have been developed by number of investigators [1,2,17]. 

Graham et al. [ 2 2 J have presented a sophisticated mechanical 
analog, for a rectangular container under lateral excitation, 
consisting of a fixed mass and infinite number of spring 
supported masses. This model has since been used extensively, 
specially in aircrafts, to study the dynamic behavior of 
fuel motion. Another mechanical model based on a fixed mass 
and infinite number of oscillating pendulum masses is reported 
in [29] for a circular cylinder undergoing lateral excitation, 
ihe advantage of the latter analogy is that the parameters 
of the model do not change with change in vertical acceleration 
Experimental verifications of the theoretical results based 
on the models have been also reported [1,2,17]. 

A mechanical model, for liquid sloshing under vertical 
or simultaneously under lateral and vert? cal excitations of a 
container, is required for vehicle response studies. Such 
a model is not available in the literature [1]. 

1 . 2 . 1 . 3 Stability of free-surface of liquid 

When a container filled with liquid, such as water, 
undergoes vertical periodic excitation, a pattern of 
stationary or standing waves can be observed on the free 
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surface. The governing equation of motion of free surface 
is characterized by parametric excitation and in a rigid 
container is described by Mathieu’s equation [8,1]. T Ae 
standing wave may have subharmonic, harmonic or superharmonic 
frequencies of oscillations. The stability of the free 
surface is accordingly influenced. 

Workers in the past have studied the stability of 
free surface under vertical excitation mainly tor rigid 
containers [1,2,17]. Benjamin and Ursell [8] have tried to 
generalize the linearized theory of stability of free surface 
under vertical periodic acceleration in rigid containers. 

For a particular case of circular cylinder, they also 
conducted experiments to confirm, the theory. It is of 
interest to extend the linearized theory of stability of 
free surface to containers with flexible walls. 

1. 2.1.4 Damning, of liQuid motions in r igio — c ont ainers . 

Damping has a profound effect on control and stability 
of liquid free surface. This is particularly true when the 
control system natural frequencies are close to the liquid 
natural frequencies. The damping of liquid, is affected by 
the liquid height, liquid kinematic viscosity, size and the 
kind of wall prepared for the container, liquids, like 
water, have negligible damping in containers with smooth 
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walls [1,8,16]. Theoretical and experimentally verified 
damping law of liquid motions are reported for containers 
with smooth, baffled and compartment ed walls [1,2,16]. 

1 . 2 . 1 . 5 Fluid-elastic interaction 

Due to external excitations, the generated hydrodynamic 
pressure deforms a container with flexible walls, which in 
turn affects the hydrodynamic pressure. This kind of fluid - 
elastic interaction occurs in several ways, such as, bending 
and breathing motion of container bottom and side walls. 

In some studies, a simple beam bending mode is assumed 
for container walls and the sloshing is treated as a forced 
motion problem [1,7,17,35]- Both theory and experiments 
confirm that the presence of free liquid surface results in 
an increase of wall bending resonant frequency over the 
container of a capped liquid. On the other hand, the flex- 
ibility of the bottom of the container reduces the natural 
frequencies of the system [l,7,ll]* 

Studies on fluid-elastic interaction have been carried 
out by assuming liquid as compressible also. Such studies 
include breathing vibration of containers [10,12,43] and 
lateral oscillations of liquids [9]- workers of reservoir 
and reservoir dam systems have also considered the fluid- 
elastic interaction problem [15,38]. 
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1.2.2 Vehicle dynamics 

9 

The trend in the design of transportation systems is 
increasingly towards achieving shortest travel time with 
maximum comfort for passengers and transport of maximum volume 
oi cargo with safety and ease. To meet these requirements 
vehicles are required to provide better riding quality, 
higher velocity and larger capacity to carry The load. For 
the design of ground based vehicles, therefore, a better 
understanding of the dynamic behavior of the vehicle rolling 
over the ground surface is of considerable interest. 

Since ground unevenness is the source of excitation 
for the rolling vehicles, the characteristics of ' the prepared 
ground surfaces such as roads, runways, railway track, etc. 
have been studied by many workers. Surface profiles are 
generally assumed to be one dimensional homogeneous Gaussian 
stochastic process in space and can be represented through 
their power spectral density. Such profiles include 
runways [23,51], highways [18,19,46], railway tracks [2l]* 

Equations governing the vehicle dynamics may have 
nonlinearities in inertia, velocity and restoring force 
terms. If the vehicle operates in random environ- 
ment, it is possible to approximate the equations of motion 
in a linear form by approximate methods. These methods 
include Fokker-Planck equation, small perturbation approach [46] 
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and equivalent statistical linearization method [14,32]. 
i’he equivalent statistical linearization method has been 
proposed and extensively used by Caughey [14] for a single 
degree of freedom model [4,26,32]- Later this method has 
been extended to multi-degree-of-freedom system by Iwan and 
Yang [26], Another method applicable to multi-degree-of- 
freedom system is proposed by Atalik and Utku [4j. Limita- 
tions of each of these methods have been discussed in detail 
by Atalik and Utku [4]- Iwan and Yang's method is suited 
for spring-mass systems, but fails to linerize inertial 
nonlinearities, ihe more direct and simple method of 
stochastic linearization is the one suggested by Atalik 
and Utku [4]. By this method a multi-degree-of-freedom 
system in inertial, velocity and restoring forces can be 
linearized. 

Analytical investigations for vehicle dynamics are 
possible both in time and frequency domain [51]. As the 
profile unevenness is assumed homogeneous, the input excita- 
tion process becomes stationary for vehicles moving with 
constant velocity [48,51]- Bor the linear model of the 
vehicle, the response becomes a stationary process. Both 
heave and pitch degrees of freedom have been considered for 
single point and two point excitations respectively by these 
investigators . 



For an accelerating vehicle rolling over an uneven 
surface, the input excitation becomes non-st at ionary. In 
Refs. [ 46 , 48 , 51 ] » the response of accelerating vehicles has 
been presented. Ref. [51] has presented a frequency domain 
approach using the concept of evolutionary spectra suggested 
by Priestley [39 ] . 

lo analyze properly the dynamic performance of the 
vehicle carrying liquid as cargo, it is desirable to consider 
the sloshing aspects of the liquid in a container [1] . The 
dynamics of such vehicles is of interacting nature and is 
difficult to carry out. However, if one replaces the sloshing- 
liquid by its equivalent mechanical model, as discussed 
earlier, then the dynamic analysis of the vehicle system can 
be simplified. Also, by selecting only a few significant 
oscillating masses of the model, the analysis is further 
simplified. One can choose an analytical method to solve 
the equations of motion for vehicles, such as a railway 
wagon carrying a liquid cargo, as discussed above. 

1 . 3 Objective and Scope of the Present Study , 

The objectives of the present study can be broadly 
classified as: 

(i) !To study the stability of the plane free surface in 

a flexible container under periodic vertical excitation i 
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(ii) to develop a mechanical model for- liquid sloshing in 
a container under the simultaneous action of vertical 
and horizontal excitations, and. 

(iii) to study the dynamics of the vehicle carrying liquid 
as cargo and comparing the results with equivalent 
solidified mass of the liquid carried by the vehicle 
under the same random environment.. 

In the first part of the work, an incompressible 

\ 

inviscid fluid in a rectangular container with flexible 
bottom is subjected to longitudinal periodic excitation, 
linearized equations of motion and boundary conditions for 
the fluid oscillations are assumed . The flexible bottom 
plate undergoes elastic deformations under thin plate 
approximations and is assumed to be simply supported, on all 
edges. The mathematical treatment for the coupled differen- 
tial equations of motion is carried out by the method 

v 

suggested by Bolotin [13]. This has bell presented in 
Chapter 1. 

In the second, part of the work, described in Chapter 2, 
a mechanical model representing the liquid sloshing behavior 
in a rigid rectangular container is presented. The container 
is subjected to simultaneous action of a constant horizontal 
acceleration and a periodic longitudinal acceleration. 
Linearized equations of motion and boundary conditions are 
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assumed for the inviscid, incompressible fluid* undergoing 
the presumed excitations. The equivalent mechanical model 
consists of a rigid mass and infinite number of oscillating 
pendulum masses. The oscillating masses are found to be 
independent of the parameters of excitations involved. 

Each oscillating mass is representative of the corresponding 
fluid mode. The fixed, mass represents the mass of fluid 
which remains unaffected by the applied excitations. The 
analogy is valid for stable free surface conditions. 

Chapteis4 and 5 present the results of the investiga- 
tion of the dynamic response of a railway wagon carrying 
liquid as a cargo. The response is determined for the wagon 
moving with constant velocity and constant acceleration. 

Chapter 4 describes the simulation of' the track 
unevenness on a digital computer and determination of the 
response of a non-linear vehicle-fluid system model through 
numerical integration. The mathematical expression, for 
power spectral density of railway track, given in Ref. [21 ] 
is modified and tested. An analytically convenient mathe- 
matical expression is chosen such that it closely approximates 
the given power-spectral-density of Ref. [21] by least-square- 
fit procedure. The random process is generated by Monte-Carlo 
method suggested in [ 44 , 45 ] using the given power spectral 
density. The power spectral density of the simulated process 
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is determined and tested. It is found that the target and 
simulated spectra are close to each other. 

ihe vehicle is assumed to respond under heave mode 
only. The liquid cargo is replaced by an equivalent 
mechanical model described in Chapter 3. The governing 
equations of motion of the vehicle system are described 
through 2-nonlinear coupled differential equations. 
Equations of motion are solved numerically by Hunge-Kutta 
algorithm. The input process is generated as described 
earlier. It is found that the vehicle response is a narrow- 
band nonst ationary random process. The coupling effect 
between the vehicle response and liquid sloshing is found 
to be negligible. The liquid response is determined for 
three different heights of the liquid in the container. 

Eor different heights of the liquid it is found that the 
relative wave height of the liquid reduces with increase 
in the liquid depth in the container. 

Chapter 5 describes the vehicle cynamics for a model 
with heave and pitch degrees of freedom. The vehicle 
dynamics is influenced by the sloshing of liquid. The 
dynamic behavior of liquid is represented through the 
mechanical analog developed in Chapter 3. The power 
spectral density of the track unevenness is assumed as in 
Chapter 4. The equations of motion describing the vehicle- 
system dynamics are nonlinear and coupled. They contain 
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inertial, velocity and restoring force nonlinearities, The 
equations of motion are linearized by the method suggested 
by Atalik and Utku [4]. ibis method, as pointed out 
earlier, has an advantage of linearizing the nonlinearities 
in inertial terras as well, ihe linearized coupled equations 
are decoupled using the method of Foss [20]. The uncoupled 
equations are transformed from time domain to space domain. 
Finally the solution is obtained in frequency domain using 
the theory of random vibration. This is possible as the 
ba.se excitation process is homogeneous in space domain. 

The response statistics for displacement and velocity are 
determined as a function of vehicle velocity. The response 
statistics are compared with the vehicle system which has 
an equivalent solidified mass of the liquid. Nine different 
vehicle- container configurations are considered for the 
study. 

The vehicle chosen is a railway wagon which is 
assumed to move with constant velocity arid constant accelera- 
tion. The base excitation to the vehicle is provided 
through the rail-road roughness whose form has been described 
in Chapter 4 • 

For constant velocity case, the response statistics 
are obtained, in analytical form. The results obtained are 
presented and compared with the equivalent solid mass- 
vehicle system. It is found that the response statistics 
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vary rapidly in the neighbourhood of small vehicle veloci- 
ties. ihe response settles down with increase in velocity. 

At large vehicle velocities the response of the vehicle, in 
heave and pitch modes, tends to be independent of vehicle 
velocity. Similar characteristics are found for the free 
surface response. 1‘he relative wave height of the fluid 
at the free surface increases with decrease in fluid depth. 

It further increases with the increase in eccentricity of the 
container c.g. placed with regard to vehicle c.g. Similar 
observations are made for the shear force and bending moment 
developed by the sloshing liquid at the bottom edge of the 
side wall of the container. 

the response of vehicle-fluid system with that of the 
vehicle-solidified mass system are compared. It is found 
that the sloshing of liquid causes reduction in the vehicle 
response over most of the. velocity range of interest. At 
some velocities, the sloshing of liquid causes increase in 
vehicle response. The range covered in such cases is narrow 
and towards small velocities. 1'he decrease in vehicle 
response is due to absorber effect, ihe liquid absorbs 
the energy from the vehicle, thus reducing the response, 
which in turn affects the liquid sloshing. Ihe absorber 
effect is found to increase with increase in eccentricity 
of the container placed with respect to the vehicle c.g. 
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For the constant acceleration case, the governing 
equations of motion for vehicle-fluid system and vehicle- 
solidified mass system are solved numerically using the 
simulation technique as described in Chapter 4 for the 
same rail-road roughness. It is found that the overall 
vehicle and liquid response is similar to the constant 
velocity case. 



CHAPTER 2 


PARAMETRIC INSTABILITIES OE THE PLANE ERIE SUREACE 
OE A LIQUID IK A CONTAINER WITH ELIGIBLE BOTTOM 
UNDER VERTICAL PERIODIC MOTION 


2.1 Introduction 

Vertical vibrations of a. container filled with 
liquid occur in several cases. In rockets and missiles, 
vertical vibrations occur due to dynamic coupling between 
the structure and the engine thrust. In aircrafts vertical 
excitation of fuel container occurs during taxi, landing 
and take-off and also during flight because of air turbu- 
lence. Containers filled with liquid in automobiles are 
also subjected to vertical excitation due to road roughness 
Other examples of interest are containers subjected to 
vertical excitation during earthquake and include dams, 
overhead tanks etc. In vertically excited containers 
patterns of standing waves are observed. At certain 
forcing frequencies the surface may respond in a large 
amplitude standing wave and sometimes may disintegrate, 
forming a dense spray. It has been also observed that at a 
large input levels of vibrations, the surface motion 
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becomes violent and bubble formation may result. Litera- 
ture on such behaviors of sloshing liquid has been widely 
reported [l, 2,3,17]. 

In a prior investigation [8], the linearized theory 
of liquid motion is developed to study the stability of 
the plane free surface in rigid containers. In the present 
analysis the linearized theory of liquid motion is extended 
to study the stability of the plane free surface in a 
container having flexible bottom and rigid side walls. 

2.2 Mathematical formulation of the Problem 

Consider a rectangular container of length, a, breadth, 
b and depth, h^ (Pig. 2.1). The container is partially 
filled with an incompressible, inviscid liquid to the height, 
h as shown in the Pig. (2.1). Dimensions of the container 
are such that surface tension is neglected. The side walls 
of the container are rigid while the bottom is flexible 
and is assumed to be simply supported on all edges. The 
container is excited vertically upwards with acceleration 
g cos ut, where is constant. The fixed coordinate 
axes (x,y,z) are attached to the bottom corner of the 
container. Within the liquid the velocity potential 
must satisfy the Laplace's equation, viz., 
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afsL + ah + Jds = o (2.1) 

o o o 

ax 3 y 3z 

i'he velocity in the three directions x,y, z are given as 




(2l 11 M) 

K 3x ? 3 y 5 3 z' 


( 2 . 2 ) 


2.2.1 Boundary conditions 


On the free surface of the liquid two conditions are 
to he satisfied. The first is the dynamic condition which 
demands that the pressure on the free surface must equal 
the pressure outside the liquid, which is assumed to he zero 
without loss of generality. Mathematically, it is given 

by [8] , 

[|f + (g - g 0 cos w tk] z=li = 0, (2.3) 

where z = C (x,y,t) defines the free surface and g is the 
acceleration due to gravity. And the second condition is 
the kinematic which requires that a particle on the free 
surface moves with the velocity of the free surface. Hence 



2L 1 

3t-* z=h 


(2.4) 


In deriving last two equations products of small quantities 
are neglected [8j. Hence the formulation given here is 
restricted to small motions. 

The normal component of velocity to the rigid wall is 


to be zero . Hence 



at x = 0, a 
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39 
3 x 



(2.5) 


3 y 
3 y 


0, 


at y = 0, 1 


( 2 . 6 ) 


Finally, if the flow does not cavitate, the particles 
of the liquid in contact with the flexible bottom and the 
bottom itself must move with the same velocity. If ¥(x,y,t) 
denotes bottom plate displacement then the last boundary 
condition is given by 
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3z-*z=0 3t 


(2.7) 


In deriving the last boundary condition the. products 
of nonlinear terms are neglected and the liquid velocity is 
evaluated at the undisturbed position of the bottom plate. 
1'he approximations are consistent with those of free surface 
conditions . 


2.2.2 Equation of motion of the bottom plate 

The differential equation for the motion of a thin 
elastic rectangular plate is 


4 4 

D [a_| + 2 _1_1 


3x 


y 2 


3 ^ 

+ * — j] + cW + m w 

3y p 


-p c ( x ,y» t ) 


( 2 . 8 ) 


where , 


¥ = transverse deflection of the plate, 
D = Et^/12( 1 - V 2 ) , 


(2.9) 
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E = Young's modulus, 

tp = plate thickness, v is Poisson's ratio, m^ is 
the mass per unit area ox the piste, c is damping constant . 
p c (x,y,t) is the total pressure acting on the plate and is 
given by 

P c U,y,0,t) =ph (g-g Q cos oj t ) - p 1^-] (2.10) 

Phe last term in this equation is referred to as dynamic 
pressure [8] and. P is the density of the liquid. 'the 
right hand member of (2.8) has a negative sign because 
positive pressure acts on the plate in a direction opposite 
to that of the positive plate deflection. Here simply 
supported bottom plate on its all edges is considered 
although the formulation holds good for other ed.ge condi- 
tions. I'or a simply supported rectangular bottom plate 
the displacement W(x,y, t) is given as follows, 


W(x,y,t) = 11 T] (x).* £ (y) q (t) 

r s 


where n r ( x ) 



^(y) = sin (~y) , 


( 2 . 11 ) 


( 2 . 12 ) 


are the normal modes in x and y directions and q (t) is the 
generalized co-ordinate for the modes r and s. In this 
co-ordinate system the differential equation for the plate 
vibrations is characterized, as follows. 
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\ s [q. rs (t) + 2 P r «. 

a b 


^ “ T.« O*)] = - p rs (2 * 13) 


where, M rg = J J sirt ~x • sin^fp y) dxdy 


P 


= m 


o o 
ba 


a 


b 


P 4 

0 = c J J sin 2 (~x) 3in 2 (~y)dxdy 


rs 


o o 

cba 
4 * 

C 


rs 


rs 


2I'i rS w rs 


(2.14) 


w rs ~ 71 LV a' T 


[(f ) 2 + 

P 

P; s (t)= J | p c (x,y,0,t) sin ~x • sin | ~y dxd.y 


S\ 2l D 


r-rc. 


sir 


2.2.3 Solution of the coupled boundary value problem 

Following functions [8] satisfy the boundary condi 
tions (2.5) and (2.6), 


cp (x,y,z,t) = 

I 

i 

l 

3 

cos (~x) cos(^y) 

d 


(2.15) 

S (x,y,t) 

I 

I 

, in % jit . ■> 

cosi— X) cos(£-j; 

j( z,t) 

(2.16) 


i 

j 



i‘he functions F. .(z,t) and Q. •( z,t) arc to be determined. 

1 J X j 

On substituting (2.15) into (2.1), (2.7) one gets 
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3 2 F. .(z,t) 

iJ- - kj. P , .(z,t) = 0 (2.17) 

9 z ^ 1 J 1 J 

where k^ = r 2 [ (J) 2 + ( £) 2 ] (2.18) 

and 

si 1 . . 

I l oos (|2x) costly) 2=Q - if(x,y,t) (2.19) 

1 (tj 


i‘o make the last boundary condition (2.19) applicable, one 
multiplies both sides of it by cos(—x) and cos( r—y) and then 
integrates with respect to x and y between the limits of 
x = (0,a) and. y = (0,b) respectively, finally one gets using 
( 2 . 11 ) 


3f ,(z,t) 

i-si - 

8 z 


z=0 


i 


y c rs 4 Q (t) 

£ ij rs v ' 


where . 


.,rs 


I6rs 


* 2 (r 2 -i 2 )(. 2 -J 2 ) 


( 2 . 20 ) 


Relation (2.20) is valid for odd-even combinations of (r,i) 
and (s,j). Por all other combinations of (r,i) and (s,j) 
one has 

P. • ( z , t ) "1 . = 0 (2.21) 

3 z ij ' ' J z=0 

i'he other boundary condition is obtained from (2.3) and (2.4) 
using (2.15) and (2.16) 
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and 


3P ,(z,t) 

U_ 1 

9t -l z=h 


+ (g - g 0 cos ut) Q ± • ( t ) = 0 


9 F ,( z,t) 

iJ ] 

9 z ■* z=h 


dt 


( 2.22) 


(2.23) 


The boundary condition (2.21) corresponds to the case in 
which there is no coupling between fluid and elastic modes 
of plate. The solution of (2.17) corresponding to boundary 
conditions (2.21) and (2.22); (2.23) has been already dealt 
with [8]. Here solution corresponding to the coupled 
boundary conditions (2.20) and (2.22), (2.23) is treated 
for (2.17) . 


The general solution of (2.17) is given by, 

F ± j(z,t) = A j_j(t) cosh(k ± jz) + B i ^(t)sinh(k ± ^z) ( 2.24) 
Boundary condition (2.20) gives 


B u (t) = 


E s c_ rs , a (t) 


L rs 


(2.25) 


1 J 

b.c. (2.22) and (2.23) determine respectively 


A. .(t) cosh (k -h) + B. .(t) sinh (k. „.h) 
id ij J- j 

+ (g-g 0 cosut) Q ± j(t) = 0 (2.26) 

and 

k ij A ij (t) 311111 (k ij h) + k ij B i.i (t) cosh (k ij h) = hj (t) 


(2.27) 
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To determine Q. .(t) one needs to eliminate A. -(t) or 

ij 

its derivatives from ( 2 . 26 ) and. (2.27). By multiplying 

( 2 . 26 ) and (2.27) with k. . sinh (k. .h) and cosh (k. .h) 

J i u i J 

respectively and then differentiating (2.27) w.r.t. t one 
gets after subtraction from these two equations 


Q ij + (g~g 0 cos w t) k ± . tanh (k^h) Q ±J (t) = 
Substituting for B. .(t) using (2.25) one gets, 

X J 


k. .B. .(t) 
x.i i.r 

cosh(k. .E.) 

J 


Q ± j(t) + k ± ^ tanh (k ± ^h) (g-g Q cos u>t) Q^t) 


l Z 
r s 


p rs 

Til 




cosh (k. -h) 
U 


( 2.28) 


Using specified initial conditions equation (2.28) is 

solvable for Q. .(t) provided q (t) is known. Further Q. .(t) 
1 J rs 1 3 

and q rg (t) are coupled. At this point, the velocity 

potential for all the plate and fluid modes is written using 
(2.15), (2.24) , (2.25) as follows. 

9 (x,y,z,t) = I l cos(~x) cos(f^y) [A. .(' t ) cosh(k. -z) 

± J ^ ^ ' 1 J 1 J 

r rs 6 ( t ) 

S. | C i j q rs (tJ 

+ v sinh ( k i i z )] (2.29) 

id J 


in which A. -(t) is known either from (2.26) or (2.27) on the 

X J 

assumption that Q^(t) is also known from. (2.28). The func- 
tion q (t) is determined using (2.10), (2.13) and (2.14). 

xo 

Using (2.10) and (2.29) cue gets, 
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' P c U,y,0,t) = p h(g-g o cos 6>t) - pi I cos(|£x) 


1 .1 


cos(^y) • iL^t) 

in which Jl ^(t) is known either using (2.26) or (2.25) and 
is given by 

Aij(t) = -[(g-g n cos ut)Q s H (t) + — 


r* r» Z*S 

S 2 r : . . q 


1 


ij 


k . 

i , 


slnh -cbSE(is,ii) 


ihus , 


P c (x,y,0,t) = +ph( g-g Q cos cat) + I cos-x-cos ^y 


^X C Z ' d r,sJ sinh k ii 


[(g-g 0 cos w t ) Q ± -j ( t ) + 11 ^ 11 


i 3 


h 


1 


k id 


1 =, 

-* cosh(k. -h) 

X J 

(2.30) 


Using (2.30) one gets from (2.14) 


Pjg(t) = P - - | ab [(-l) r -l][(.-l) S -l](g-g cos u t) 
■re rs 


+p | ( r k~^Sitb~~hT [ k ii ( e-'So GOS “*) Q ij (t) 


C rs 


+ siirn (k^fc) ^ o. k 1 - 


q r s 
r l s l 


(2.31) 


On substituting (2.31) into (2.13) one gets 
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14 + 28 03 4 (t ) + 03 2 g (t)] 


ph ab ((-l) r -l)((-l) S -l) 

= - [ 5 (g-g rt cos tut) 


2 

ix rs 


+ D y y ab ~xj 

f 4 k. . cosh(k. .h) 
1 J xj ID 


C rs 


X* s 

4 sinh k. -h £ £ c.^ ^ q 

lj r 1 s 1 ij r i s i 


(\j(g-g 0 co s o)t)Q i ^(t) 

(t))] (2.32) 


Rewriting for convenience, Lq. (2.28) one has 


H C iD a rs (t) 

Jl o 


e ij (t) + k ij tsmh ( k ij h )-(e-So cos = - Tolh '- (k ~T 

X. J 


(2.33) 


Last two equations are sufficient to determine as 

many unknowns as in Q. i (t) and q (t) depending upon the 

initial conditions specified for the same number of unknowns. 

1‘hus, the problem is well defined for the unknowns Q. .(t) 

x D 

ana q rs (t). It has been also noted earlier that (i,r) and 
(Djs) individually should have odd-even combinations so that 

ICS 

0. . is other than zero for a coupled boundary value problem, 
x D , 

X* s 

For other combinations, C. . does not exist and belongs to 

-L ti 

uncoupled, boundary value problem. It is clear that the 
solution involving axl the fluid and elastic modes is 
difficult if not impossible, therefore to gain insight, 
the problem is treated for one elastic and one fluid mode only. 
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One of the many possible combinations of (i,r), (j,s) may 
be either r=s=l for plate and i=j=2 for fluid mod es , or, 
t=s=2 for plate and i=j=l for fluid modes. 


Presently as an example, r=s=l and i=j=. 2-are assumed. 
Further for simplification the plate damping, P rS » is assumed 
zero and hence one gets from (2.32) and (2.33) the two 
coupled, equations as follows. 


fc%i +— TT 


P ab( C 22 ) 2 tanh k 22 h 


22 


) ( "k ) I'-i 1 w 1 1 Q.-} 1 ( ) 


ii w ii h h' 


/ 4 phab pab °22 

- ( ~ 2 + 4" cosh (k^h) 1 ( «-4o cos 


71 


(t) 

Q„(t) + k 22 tanh(k 22 h).(g-g o cos cot)Q 22 (t) - 


On rearranging the last two equations one gets 

• « A 

9-11 + “pi Q.pp (iO + D 2 (l-e cos oj t)Q 22 (t ) = -D-^l-e cos at) 

(2.34) 

“ D 3^11^ + '^22^ +n 22^ 1_£ C0S Q 2 2^ = ° (2.35) 


where, 

* 2 
“ 11 


M. 


M. 


Cl 22 = g k 22 tanh (k 22 h) 


11 


M 11 + 


p (cio) 2 *ab tanh( k 00 h) 




22 


4 k 


22 


41 _ 16_ 
' 22 = Six' 


e = 


g 


(2.36) 
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D- 


D, 


4 hpab 

2 * 

* M 11 
r ll 

_22 

cosh( k" 22 h) 


g; -D o — 


P c 22 ab 

4M£-^ cosh( k- 22 h) ® ’ 


(2.36) 


Combining equations (2.34) and (2.35) in the matrix form 
as follows. 


C I' + (l-A-sB cos lit) P = D 
where 

C ~ r 1 t r 

— ~~ J » I - [q l] » 


(2.37) 



*2 

*1 

f 


— 

—1 

A = 

i -“n 

0 

- E 2 

1-^22 1 

I 

1 9 

| 

B = 

j 

0 

0 

2 

a 22 


— 

J 

1 


— 



(2.38) 


P = i 


H 11 i 

^22 f 


D = 


-D-j-( 1-e cos w t ) j 
0 


Equation (2.37) will be used to study the parametric inst- 
abilities of the liquid surface. 


2 . 3 The Parametric Instabilities 

The parametric instabilities of a liquid surface 

mainly differ from the forced vibration response in 

\ 

following respects. 

(i) The parametric instability may occur at frequencies 
other than the natural frequency of the free-surface oscil- 


lations. 
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(ii) For transverse excitation of tank, the frequency 
of the free surface is precisely the same as that of the 
forcing frequency and resonance occurs at discrete frequen- 
cies. In contrast to this, the parametric resonance exists 
over a continuous regions of excitations (regions of dynamic 
instabilities) in which we are interested [13] * 


Parametric instabilities are classified (ljj as primary 
and secondary instabilities. If a system has natural fre- 
quencies ft n (n = 1,2,3,...), then the primary instabilities 
2 

occur at w = wi^k & = 1,3,5,..., as the excitation 

parameter tends to zero, where w is the frequency of periodi- 
cally varying parameter. Similarly, the secondary instabi- 

2 

lities occur at w = r- ft with k = 2,4,6,..., as the excitation 
parameter tends to zero. Further distinction in parametric 
instability is made when k=l and 2 and corresponding inst- 
abilities are termed as principal primary and principal 
secondary instabilities respectively. Bolotin [13 ] has 
suggested that principal instability is the most critical 
one among all others and has greatest physical importance. 
Principal primary instability region is also known as the 
principal region of dynamic instability. 


2.3.1 Determination of regions of dynamic instabilities 

io study the dynamic instability the homogeneous part 
of equation (2.37) is considered viz., 
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CF + (I - A - eBcos u t) P = 0 (2-39) 

Bolotin |3L3] has shown that the regions of unbounded 
solutions are separated from the regions of bounded, solutions 
by the periodic solutions with period i sue 21, where 

O ■rr 

i = — (Ref. Appendix A). More explicitly, the region of 
instability is bounded by two solutions of identical periods 
whereas the region of stability is bounded by two solutions 
of different periods, ihe regions of primary instability 
are contained by the solutions with period 2i and secondary 
instability regions are contained by the solutions with 
period 1. 

2 . 3 . 1 . 1 Primary instability regions 

To determine the regions of primary instability of 
the system (2.39) one assumes [13] 

P = l ^ [a k sin (| u>t) + b k cos (|ut)], (2.40) 

k=l ,3,5 

where a k and b k are unknown vectors and are independent of 
time. Substituting (2.40) into (2.39) and comparing co- 
efficients of sin(^~-) and cos(~^) on both sides of the 
equation one obtains the following system of matrix equa- 


tions. 
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(I - 

A 

+ 1 

e B - 

05 

4 

C) 

*1 - 

1 

2 e B a^ = 

= o, 


(I - 

A 

1 

* 4 

k 2 m 2 

c) 

£k 

1 

J_ 

2 6 

- k-2 + 

— k+2^ 

= o» 

(I - 

A 

1 

- 2 

s B - 

1 

4 

05 2 ( 

D 

- 1 e ® 

= o, 


(I - 

A 

1 

” 4 

, 2 2 
k 05 

c) 

*k 

1 

2 £ 

® <k_ 2 + 

-k+2^ 

= 2* 

with 

k 

= 3, 

5 f • » < 






(2.41) 


Combining all these infinite set of matrix equation one gets, 



etc. 

The system of equations (2.42) has solutions other than 
trivial provided the determinant of the co-efficient matrix, 
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[ G-] , is zero. This will ensure the determination of frequen- 
cies bounding the instability regions, The determinant, |G] , 
is of infinite order. It can be shown that, |&|, belongs to 
the class of normal determinant and is absolutely conver- 
gent [13]. i'he direct application of j&j to determine the 
exact instability regions is practically inconvenient. 

Hence, approximately the regions of instability are obtained 
by truncating the series in equation (2.40) at k=l . This is 
equivalent to consider the determinant circumscribed with 
dashed lines in (2.42). The approximation, k=l, necessarily 
yields the principal instability regions. For a better 
approximation and also to determine the higher order regions 
of instability one has to truncate the series in equation 
(2.40) at higher values of k. Bolotin [13] has shown through 
an example of a beam subject to a periodic axial force that 


t The 

determinant 





1+C 11 

j- 

C 12 

°13 

0*0 

A = 

C 21 

1+0 22 

c 23 
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1 + C-Z >7 

33 

O * * 

is called normal 

if the 

double series 


abs o lut ely convergent . 


I 1 

i=l k=l 


°ik is 
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for values of e upto 0.6, the regions of principal primary 
instability obtained by k=l approximation, are within one 
percent of those obtained by k=3 approximation. In another 
example of coupled equations Bolotin has shown that the 
width of the regions of dynamic instability rapidly decreases 
as the number of regions increases. The width of the inst- 
ability region reduces as s r , where r is the number corres- 
ponding to the first, second etc. instability regions. Thus, 
the first region of instability has the greatest width. 


2 . 3 . 1 . 2 Secondary instability regions 

To determine the regions of secondary instability, 
which are separated by solution of period i, one assumes [4-3] 

2 = I h> + k | ? _ 6 ts-k sin + it 003 (2 - 43) 


Substituting this equation into (2.39) and comparing the 

yields the equation 


co-efficients of sin(lpt) and cos(~t) 
of the form (2.42) with the column vector 


a, 
i bi 
a' 


matrix [G] having the elements as follows. 


and the 




G H = 1 ~ & 12 = °» G i3 = -eB, & 14 = 0 

G 21 = 0, & 22 = l ~ A - w 2 C, G 23 = 0, G 24 = - | s B, 

= 2 " a B» G'J2 =: ^ , ^33 = — ~ 6,* G ^ 4 =l J * 

^ 4 ]_ = ^42 = — 5 G 43 = ^ ’ ^44 = ( L ~ A - 4w^C), ..... 

et c . 
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The determinant corresponding to k=2 approximation is 


'll 

G 12 

G 13 



21 

G 22 

G 23 

■ 

(2.44) 

'31 

G 32 

G 33 

i 

f 

! » 



and when equated to zero, determines the regions of secondary 
instability approximat ely. 

2-4 Results and Discussions 
2-4-1 Computations performed 

Using the analysis presented in the earlier sections, 
the following computations have been performed. Regions of 
primary and secondary instabilities associated with the first 
mode of the plate (1,1) and second mode of the fluid (2,2) 
are determined so that coupling effect is retained. The 
instability regions have been presented as plots of the non- 
dimensional frequency versus the excitation parameter e. 

The non-dimensional frequency is the ratio of exciting 
frequency to fluid natural frequency associated with the 
fluid mode under discussion. The excitation parameter is the 
ratio of vertical acceleration to acceleration due to gravity. 
Values of e up to 0.9 have been considered for both primary 
and secondary instability regions. 

The regions of instability have been obtained for three 
ratios of plate to fluid natural frequencies. These ratios 
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have been chosen as to cover the range of interest 0.965, 

1.07, 0.512. 

For the determination of primary instability region 
equation (2.42) with the circumscribed matrix and for 
secondary instability region expression (2.44), by equating 
it go zero, are assumed. 

2*4.2 Effect of flexibility of container wall on stability 

As expected the effect of flexibility of the container 
bottom is to narrow down the region of instability with 
respect to the rigid container. This effect is markedly seen 
on the principal secondary region specially towards the higher 
values of excitation parameter e (Fig. 2.2) . At the lower 
value of e the effect on instability region is marginal as 
compared to rigid container instability region. 

On the principal instability the effect of flexible 
bottom is to introduce one more instability region as 
compared to the rigid case. One of them lies near the 
plate natural frequency and the other near the fluid natural- 
frequency. This indicates the coupling action between the 
fluid modes and wall modes. Same effect is also observed 
for the ^ subharmonic (Fig. 2.2a). It is seen that for the 
~ subharmonic case the introduction of second region is 
almost invisible. It is to be noted that the introduction 
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of second region of instability i s visible only when the 
natural frequencies of fluid and plate are close to each 
other (Fig. 2. 2(a), (b)) whereas for differing values of 
natural frequencies of liquid and plate, such as shown in 
figo 2.2(c), the likely second region of instability is 
almost indistinguishable. The ratio of plate to liquid 
natural frequency is chosen for Fig, (2.2(c)) is 0.512. In 
tins particular case the principal primary region is not 
appreciably affected by the flexible bottom. whereas the 
effect on principal secondary region is seen to reduce it 
markedly towards the higher values of e . 

io sum up the flexible container wail in general 
reduces the instability region of the fluid free surface as 
compared to the rigid container. Closer the natural frequency 
of fluid surface to the container wall natural frequency, 
the greater is the coupling effect between the modes of the 
two. Correspondingly greater is the effect on the stability 
of fluid surface and the tendency is to narrow down the 
instability region as compared to rigid container. For a 
wider difference between the natural frequencies, the fluid 
suiface behaves as in a x*igi<3 containex* • 



CHAPTER 3 


A MECHANICAL MODEL E'OR LIQUID SLOSHING IN A CONTAINER 
UNDER THE SIMULTANEOUS ACTION OF LONGITUDINAL 
AND HORIZONTAL EXCITATIONS 


3 * 1 Introduction 

The motion of liquid in partially filled container 
affects the dynamics of many physical systems, such as, 
aircrafts, missiles, railway wagons carrying fluid cargo, 
overhead tanks and dams during earthquake excitation. In 
such systems., there exists an interaction between the overall 
structural dynamics and sloshing of fluid. To solve the 
complete dynamics problem, it is convenient to represent 
the sloshing liquid by a simple mechanical system. 

Some investigations have been carried out in the 
past to represent the sloshing of liquid by a mechanical 
model under lateral excitation [1,22,24,29]. Graham et al. 
[22] have given a comprehensive analogy for lateral sloshing 
of liquid by moving spring supported masses and a rigid, 
mass in a container. In Ref. [29], the sloshing liquid 
is represented by moving pendulum masses and a rigid mass 
for the lateral excitation. The choice of a particular 
representation depends on the physical problem and analytical 


convenience. 
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Ground based vehicles are subjected mainly to two 
types of excitations, i'he first is the vertical excitation 
due to surface unevenness and the second is the horizontal 
acceleration due to varying velocity of the vehicle. I'he 
dynamic analysis of such vehicles is considered mainly to 
understand the design aspects of the vehicle structure, 
passenger comfort, cargo safety and controllability of the 
vehicle. Due to travel over an uneven surface, the vehicle 
dynamics is random in nature. I'he problem of dynamic 
analysis is further complicated, if the cargo is liquid. 

I 1 his is due to interaction between the sloshing liquid and 
the vehicle structure, i'he analysis can be simplified, if 
the sloshing fluid is replaced by equivalent mechanical 
model. 

Developing a mechanical model for a generalized 
vertical and horizontal excitation of a sloshing fluid, is 
relatively complex and difficult. Under vertical excita- 
tion, the governing differential equations for the liquid 
motion are characterized by parametric excitation. I'he 
response of the liquid motion may have stable or unstable . 
subharmonic, harmonic or superharmonic nature [1,8]. 
i'herefore an equivalent mechanical model should attempt to 
reproduce these responses. A comprehensive analysis is not 
reported in uhe literature even for simple cases. An attempt 
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has been made in the following analysis to represent a 
sloshing liquid in. a rigid rectangular tank, under simul- 
taneous action of vertical and horizontal excitations, by 
an equivalent mechanical system. I* he analysis is valid for 
stable system. 


3 . 2 Mathematical formulation of the Problem 

Consider a rigid rectangular container partially 
filled, with incompressible, inviscid liquid to the height, 
h as shown in the Pig. 2.1. Dimensions of the tanks are 
such that the surface tension is neglected. The coordinates 
are moving, as shown in the Pig. 2.1, with the container 
such that the mean fluid free surface has the equation. 
z=h. The container is excited vertically upwards with an 
acceleration g Q coswt, where g Q is constant. The container 
also moves horizontally with a constant acceleration C Q . 

The Eulerian equations of motion [31] for the fluid are 
given by 
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r I‘he equation of continuity 

3 V 3 V 3 V 

3x + 3y + 3 z : 


is 


0 . 


(3.2) 


If the motion was originally started from rest, there 

is a velocity potential such that (V , V , Y ) = {-3- iL_ iL_Vo 

x ‘y z' v 3x’ ay 5 az' T 

ihe equations of motion then have the integral 


It + 2 + V y + V z } + P + (g-g 0 cos w t)lz-h)-C 0 x=P(t),(3.3) 

where F(t) is independent of (x,y,z) and may he put equal 
to zero [8,31,47]. Equation (3.2) may be written as Laplace 
equation 


a 2 ? 3^9 a 2 cp 

3x 2 3 y 2 + 3 z 2 


0 , 


(3.4) 


3.2.1 Boundary conditions 

If the equation of free surface is z = 5 (x,y f t) the 
kinematic surface condition is 
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If the deflexion and slope of the 
then the products of V x » Y , C 
he neglected. ihus the linearized 
are obtained from ( 3 . 3 ) and ( 3 . 5 ) 


0 . 

free 


and 
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(3.5) 

surface is small, 
their squares can 
surface conditions 
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§T + (g ; g 0 cos ut)c - C Q x = 0 

3S. _ 3_9 

9 t 3 z 


(3.6) 

(3.7) 


Other boundary conditions require that the relative velocity 
no rmal to the wall be zero. Hence 


9 9 
9x 


9jP 

3y 


99 
9 z 


= 0 , at x = 0 ,a 

= 0 , at y = 0 ,b 

= 0 , at z = 0 . 


(3.8) 

(3.9) 

(3.10) 


3.2.2 Solution of the problem 


It can be proved that 5 and 9 can be expressed in 
terms of the same set of eigenfunctions [ 8 j . ihus one 
assumes that 


CO 00 


s (x,y,t) = l l r (t) 3 (x,y) 

i =0 j =0 1 J 1 J 


(3.11) 


00 00 


cosh (k. z) 

9(x,y, z ,t) = I I f. .(t) S id (x,y) (3-12) 
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It is to be noted that the boundary conditions (3.8) to 
(3.10) are satisfied. ihe free surface kinematic condition 
(3.7) is also satisfied. I'or the second free surface 
condition ( 3 . 6 ) one finds that 


11 r. ,(t) 


cosh (k. .z) 


1 ] J 

ij' w k.. ^sinh"( k!^ n .h) ®ij( x, ^ r ) + (S“S 0 ccsti) ^) r j_j("k) 


ij 


ij 


S id (3E » y) = C o x 


As the right hand side is independent of y, modes corres- 
ponding to j=0 are only possible. On multiplying both sides 


by cos(~— x) an d integrating over x= 0 ,a one finds that 

r ± (t) + gk. tanh(k.h)(l - e cosut)r. = -A. n 2 , for ± = 1 , 3 , 5 . 


40 

where, A. = — 

a • gk? 


= 0 , for i = 0 , 2 , 4 , 6 ,.. 


s ? = g k 2 tanh(k i h) 


e 



£ | < 1 


(3.14) 


i'he summation sign over j drops out. If r..(t) is determ in ed 

X 

from (3.14) with the prescribed initial conditions then all 
of the boundary conditions are satisfied, ihe resulting 
expression for velocity potential for zero initial conditions 
of r^(t) is given by 
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9 (x, 2 ,t) = I r.(t) oo s( k.x) (3.15) 
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where r jL (t) is given (see Appendix B) for small e as 
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(3.16) 


3.2.3 Determination of forces and moments 

i'he total pressure acting at any point is the sum of 
static pressure and dynamic pressure, xhis is given by 
Eq.. (3.3). After neglecting the second order terms, 3S done 
to arrive at (3.6), and putting F(t) = 0 one gets 


P = p ( g-g 0 cosmt) ( h-z) -P + C Q x p 


(3.17) 
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loual force acting in x -direction per unit width of the 
container is given by 


/ pdA 
A 


where A is the area of the wall. 


Inus 


F = 

x 


h C 


+ 2P 


i=lj 3,5 , 


r ± (t)/k£ 


(3.18) 


inhere K Q = pah is the total mass 
width of the container. 


o± the fluid per unit 


Moment of the continuum model 


The moment exerted by the sloshing fluid is determined 
about the centre of gravity (c.g.) ox the fluid. The 
counterclockwise moment is assumed to be positive when 
viewed in positive y-direction. The total moment, due to 
forces acting on the walls of the container having unit 
width, is given as 
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MECHANICAL MODEL Wl I li uiNt r lAt-u 
MASS AND INFINITE NUMBER OF SLOSHING 
MASSES REPRESENTING FLUID MOTION 
UNDER SIMULTANEOUS ACTION OF VERTICAL 
AND HORIZONTAL EXCITATION OF CONTAINEF 

MOTION • 
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3.5 Mechanical Flo del 

Ihe criteria for constructing the model to be equi- 
valent to the fluid motion in the container are chosen as 
follows. I he mechanical system must exert the same force 
and to cal. moment on the container as the fluid does for the 
type of container motion under consideration. l‘he fluid may 
have infinite number of modes . Hence corresponding to 
each sloshing fluid, mode there must be a oscillating mass. 

In this way infinite number of oscillating mechanical 
masses are required. However, a detailed analysis would 
show that the size of the oscillating mass decreases 
rapidly with increasing mode number [22,29]. 

fig. 3.1 shows a typical mechanical model represent- 
ing the fluid motion under vertical acceleration, g Q coso>t , 
and horizontal acceleration C Q . 1‘ he oscillating mass, m^ , 

is attached through a pendulum rod with length 1^„ The 
fixed mass, is placed at a distance, 1 , from the fluid 

c.g. The pendulum mass, mu, is confined to move in vertical 
plane. The pivot of the pendulum is placed at a distance, 

1. from the fluid c.g. 
i 

3.3.1 Determination of model parameters 

The length of the pendulum, 1^ , is determined on the 
basis of the assumption that the natural frequency of the 
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pendulum is t ae same a.s that of the corresponding fluid 
mode. Hence 



(3.20) 


The equation of motion of the pendulum is written as 


X i m i h + WS (1 - £ cos & t ) z ± = -n 1 l ± C Q * (3.21) 

lh.e solution of equation (3*21) can be obtained for zero 
initial conditions using the results of Appendix B. Simi- 
larity between Eq. (3.14) and (3.21) may be noted. 

The mechanical masses are determined by equating the 
total horizontal force exerted by the me chani cal model and 
continuum model. The total horizontal force of the 
mechanical model is given by 


E, 


x 


= in 


ov 


C 4 - 


I 

i=l 


m. C 
l o 


+ 


l 

i=l 


m. x. 
r i 


\(t) 


( 3 . 22 ) 


As the horizontal force for the two systems is same hence 
one gets after comparing (3-18) and (3.22) the following 


M G = ni G + I m. C 
oo ov o 10 


or, m 


ov 


M - I m. 
o . 1 


and, 


r ± (u) 


l “iU 2 i = 2p l 

i=l } 3,5... kf 


(3.23) 
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As each fluid mode is independent of' the other and corres- 
ponding to each fluid mode there is a mechanical pendulum 
hence 

r, (t) 


m. 1. 2. = 2P 1 

1 l x 


k 2 

1 


( 3 . 24 ) 


On substituting the expressions of 2. and r.(t) using 


the results of Appendix B one can have 


SI? 


“it ] C 0 = 2P -i A ± [ ] 


( 3 . 25 ) 


1'he terms in the brackets on the two sides of the equation 
are same. Hence on substituting for A ± , & ± from (3.14) 
one gets 


m i 8 tanh (k^h) 


lvi 


tt 3 i 3 R 


where, R = aspect ratio of the container = — 


( 3.26 ) 


Moment of the Mechanical Model 

The distance of the masses a.re measured positive 
a oove the c.g. Thus the positive counter clockwise moment 
about the c.g. is given by 


M, 


Ty 


[m L C + l m 1 C + £ m. 1.1.2. ] 


110 


111 


( 3.27 ) 
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If the two systems generate the same moment then on 
comparing (3.19) and (3.27) one obtains. 


m. 1. L. z. 
1111 


% B. r. 

.2 ii 


(3.28) 


vnere 


h - (1 - '- 1 ) 1 ) i %nL ( V 


k^tanhOLh 


(3.29) 


Using (3.24) ana (3.28) one 


o bx axns 


(3.30) 


Finally comparing (3.19) and (3.27) one obtains the 
fixed mass distance, L„, as 


m l + E in. 1. = j p"^ — 1 

ov o 11 12 


1 n ^ 

or, L = -=- [M 2_ 
o m Qv L o 12h 


i=l, 3,5 ... . 


nuL,] 


(3.31 


The moment of inertia, I , of the fixed mass, m 

ov 5 ov 

is determined by pitching the tank about the axis? y passing 
through c.g. of the fluid. The analysis is available in 
Bef . [22 j . Thus one obtains, I from xhe following. 


I -~ Tr + m if’ + y m. ( L . —1 . ) 2 = I r 1 — — 
ov ov 0 iii^,^.. 1 . 1 sy “ 1+R 2 

7Aft « tanh ( 2n+l) §£ 

+ L2° y _ £_ -j 

B(1+B 2 )it 5 n=o (2n+l) 5 


(3.32) 
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where I is the moment of inertia of fh , 

J inertia of the solidixied liquid 

mass about thp v ov-i « 

y axis passing through’ c.g. of the fluid. 

J-Iius all -Che parameters of the mechari r -i -,„e 

na - al Cc;.l mo d el are known o 



CHAPTER 4 


DIGITAL SIMULATION OF THE DYNAMIC RESPONSE OE LIQUID 
IN A CONTAINER MOUNTED ON A VEHICLE MOVING 
ON AN UNEVEN SURFACE 


4 .1 Int ro duct ion 

Liquid is often carried by vehicles as a cargo or fuel. 
The unrestrained free surface of the liquid in a container 
has a tendency to respond quickly even to a small motion of 
the container. Hence, the response of liquid in a container 
mounted on a moving vehicle is of considerable interest. 
Vehicles operating on the ground are subjected to base 
excitation due to ground unevenness. The undulations in 
prepared ground surfaces, such as, roads, railway tracks, 
runways are random in nature and, therefore, a vehicle moving 
on such a surface is subjected to random excitation. A 
random vibration analysis of vehicle-fluid system is therefore 
required. Usually the mathematical equations characterizing 
such a system, are nonlinear and involve time-dependent- 
parametric excitations. Due to such complexities, an 
analytical treatment of random vibration analysis is not 
always possible. In such cases, the digital simulation 
offers a convenient approach to determine the response. 
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Basically, digital simulation techniques [40 j involve 
generation of random processes on digital computer to match 
some of the statistical properties of the natural phenomena. 
The simulated, process is treated as an input to the system 
and che response is calculated by numerical integration of 
the governing equations. An important advantage of simula- 
tion tecnnique is that it can handle fairly complicated 
systems. The disadvantage lies in the fact that it may 
require large computer time. However, this practical dif- 
ficulty is overcome by present day high speed digital 
computers . 

In the following analysis, the Monte-Carlo technique 
has been employed [45 j to determine the response. The input 
base excitation is simulated corresponding to a given rail- 
road power spectral density [21], and fed to the nonlinear 
system. The response is determined by numerical integration 
of governing equations. The system is described later in 
the present chapter. 

4 . 2 Rep resentation of Bower-Spectral -Density of Rail-Road- 
Roughness 

In Pig. 6 of Ref. [21] the one sided power spectral 
density of rail road roughness has been given. The same 
is reproduced on a different scale in Pig. 4.1 for reference 




P.S.D. curve ref.gtj 
Simulated P.S.D. 

Least square fit curve 


5 one sided p.s.d.,ft^/(cycte/ft) 


C s 1.5172684 x 10 
/= 5,3633526 x 10 


FIG. 4.1 TARGET AND SIMULATE 
RAIL ROAD PROFILE. 
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Ref. [21 ] has discussed 'various forms of rail-road -spectrum. 
Following,, form of one-sided power-spectral-density is used 
for the analysis employed in Ref. [21]. 

. , , 2% A o 2 

® YY ^ ~ 7~2 TT 2 s /( ra ^/ /f ' t ) (4.1) 

(ft + ft c ) o 

where ft c = 0.2513 rad/ft, is the cut-off spatial frequency. 

is the constant which may he determined from the plot, 
i'he value of the constant, A , changes from track to track. 

'I he track roughness changes as it wears out with time and 
usage. v/hile the form of power spectral density (hereafter 
written as p.s.d.) remains the same as in (4.1), the cons- 
tant, A , increases with time. The choice of analytical 
expression to represent measured p.s.d. is not unique. A 
different expression may be chosen to suit the analytical 
convenience. In the present analysis a least-square-fit 
procedure is used to fit the following analytical form to 
the data in Ref. [21]. 

$ (ft) - — — , one sided p.s.d., ft /(rad /ft), 

1 X Y + ft 

where q in rad /ft and, 

C = (9 .533277) 10" 6 , 

y = (3.3698937) 10~ 2 (4.2) 

are constants. It is seen from Fig. 4-1 that the given 
p.s.d. curve of Ref. [21] and the least-square-fit curve (4.2) 
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are quite close to each other. It differs marginally 
towards the higher frequency range. Hence the p.s.d. 
represented by (4.2) is accepted to be used, in the subse- 
quent analysis. 


4 • 5 'target and Realized Power Spectrum 

the specified power spectral density is called the 
target power spectral density. the spectral density 
estimated from the simulated, data is called realized 
power spectral density. In the present analysis, the 
target p.s.d. is specified by relation ( 4 . 2 ). the spatial 
function Y(s) is simulated as follows [ 44 , 45 ], 


T(s) 


H 

= V'2 l 

kil 


A^ cos 


(a k s - r t ), 


(4.3) 


where 



kA Q 

(k-l)Afl 


$ yy(fi ) dQ— $ yy(^ , 


(4.4) 


and r^ are random phase distributed uniformly between u 
and 2 %. «yy(Jl) is defined by (4.2). the degree of 
approximation of the realized function (4-3) depends upon 
the cut-off frequency, Q , and whether Aft is small enough 
to hold (4.4). Generally the cut-off frequency, ft c , is 
chosen such that contribution to the total power beyond 
is negligible, the degree of approximation between the 
theoretical spatial function and. the realized spatial func- 
tion is ensured, further as follows. A series of digitized 
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sample function is generated using (4.3). ihe realized 
power spectrum is estimated from the realized sample function 
series [5,6,27]. In the present analysis the Blackman-1 ukey 
procedure of power spectrum determination is employed. r i*he 
smoothing of power spectrum is carried out using Parzen 
spectral window, ihe realized power spectrum and the target 
power spectrum are compared as in Pig. 4.1. It is seen that 
the target and. realized p.s.d* are close to each other. 

4 -4 .Response Liquid-Filled-Container .Mounted on 
Moving-Vehicle 

A rigid rectangular container filled, with liquid is 

mounted on a vehicle-wit hout-tyre which moves over a surface 

as shown in the Pig. 4.2. 1’he mass of the empty vehicle is 

m . It is supported, by a linear spring and a viscous damper 

whose constants are k^ and respectively. Ihe vehicle is 

subjected to vertical displacement Y(s) as shown in the 

* 

figure, ihe mass of empty rigid -container is assumed to be 
small as compared to the mass of the loaded vehicle and is 
neglected. ihe liquid in the container is idealized by a 
pendulum with mass m , length 1^, and a fixed mass, m Qv . 
ihe expressions of these masses and their respective posi- 
tions, as shown in the Fig. 4.2, are given by Eqs. (3.20) 
to (3.32). ihe pendulum is attached to a damper with 
damping constant, c v [16], as shown in the figure. 
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4-4-1 Equations of motion 

1'he equations of motion, for the moving vehicle model 
as shown in Eig. 4-2, are given by 


(m + m + m ) + jr 1 {'&, <=-i r , 

v e ov v' 1 v v v z 3 sir 2 


- t \ 

^ * r 2 ^ COS Z 


+ c i ( 2 i - y) + ^2. ( z i - tO 


(4-5) 


n V 1 v z 5 + “V'-v z i sin z 5 + c v z 5 + Vv e sin z 3 = u - (4-6) 


These equations are nonlinear and are coupled - 


4-4-2 Solution of the equations of motion 

Equations (4-5) and (4-6) are solved, numerically as 
follows. The input excitation, y(t), is the surface rough- 
ness of the track over which the vehicle moves. Eor the 
present analysis, the surface roughness is chosen according 
to the p.s.d. of the rail-road ro ugliness given by (4.2). 
Y(s) is generated with the help of Bq. (4.3) and the proce- 
dure outlined in Section 4.3. 'The vertical velocity, y(t), 
is given by 



where v is the velocity of the vehicle and the spatial 
derivative Y\ is known from (4-3)- With the known initial 
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conditions for 2-^, an< 3 z ^> one can solve the equations 
(4.5) and (4.6; numerically using the fourth order Runge- 
Kutta integration method. Since the input excitation is 
random in nature, it is necessary to calculate the response 
as an ensemble average of ij samples. Further averaging is 
necessary in terms of moving average over n consecutive 
step readings [6,27]. The step length corresponding to a 
suitable step length of integration may be chosen for Bunge- 
Butt a algorithm. 

4 . 5 Results and Discussions 
4.5.1 Simulation 

In Section 4.5, the procedure for simulating the 
random, process is outlined. The spacing AA is chosen for 
Eq. (4.5) in view of the approximation (4.4). There are 
two ways to determine the number of spacings H in the fre- 
quency domain. The first is to assume the spacing, A A as 
constant. The second is to assume the variable spacing. 

In the second, finer spacings are chosen in those domains 
of frequency where the fluctuation of the spectral density 
is more rapid and coarser spacings elsewhere [44 J. Number 
of simulation terms, li, will be decisively less with the 
second method. To save the computer time, the advantage 
lies with the second method. In the present analysis, the 



second method is employed to select the terms N. 

From Fig. 4.1, it is considered sufficient to choose 
the cut-off spatial frequency to be 1/8 cycle/ft. Beyond, 
this value, the p.s.d. is close to zero. Accordingly, total 
number of 51 variable spacings are considered sufficient, 
i'he random phase angles, r^, are generated by using the 
built in subroutine BAN of DSC system 1090. 

4 • 5 • 2 Pa rameters of vehicle, mechanical model and initial 

conditions 

f he vehicle under consideration is the one of Ref. [42] 
Important data of the wagon are listed in fable 4.1. For the 
sloshing of liquid, in a. rectangular container representing 
a wagon, a mechanical model is assumed. I'he mechanical 
model parameters are determined through the Eqns. (3.20) to 
(3.32). As the fundamental mode of sloshing liquid is more 
important than other modes [22], the pendulum, mass corres- 
ponding to fundamental mode is considered. i'he parameters of 
the mechanical model necessary for the present analysis are- 
given in fable 4-2. ihe damping ratio for the pendulum is 
assumed to be 0.0034 as suggested in Rex. [16]. 

Ihe vehicle is assumed to maintain constant accelera- 
tion. Results are presented for the case in which vehicle 
attains a. velocity of 150 km/hour over 1 km distance from 
rest, f he distance and velocity at any time, t, can be 
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calculated as follows. 

8 = V + | V 2 • 

v 2 = v 0 2 + 2C 0 s , 

where C Q is the constant acceleration s and v are the 
distance covered and velocity at time t. v Q is the initial 
velocity which is zero at t=0. 


(4.8) 

(4.9) 


Initial conditions 

Calculations are carried out for the zero initial 
conditions of the variables involved in Sqn. (4.5) and (4.6). 

4.5*3 Numerical integration and averaging 

The equations of motion (4.5), (4.6) are solved using 
the fourth order Bunge-Kutta, integration method. The step 
length of integration is usually taken as [6,27] 

A t = ~ , (4.10) 

where f is the cut-off frequency in cycle/sec. In the 
present case, the cut-off spatial frequency considered is 
1/8 cycle/ft. Hence the step length of integration in 
spatial domain is As = 4 ft . Hence 

At = — second, (4.11) 

where v is the velocity in ft /second. In the present analysis 
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the vehicle has variable velocity run. Hence a suitable step 
size is required to integrate the equations of motion. 

Prom fable 4.2, one notes the following. For all the 
heights of fluio. considered, m the container, the completely 
full container offers the least ratio of pendulum to vehicle 
natural period and is found to be approximately 5.7566282. 
Accordingly, the step size, At, equivalent to one-tenth of 
the vehicle natural period, is considered sufficient in the 
present analysis. 

Since the input excitation function is random in 
nature, it is necessary to take the ensemble average of 
sample functions. It was found sufficient to take ensemble 
average over five sample functions, further averaging is 
carried out over the five consecutive readings of the sample 
response function. 

4.5*4 Vehicle response 

Vehicle response with varying fluid height in the 
container is studied, ihe mean squar response of vertical 
displacement with non-dimensional time is presented in 
Fig. 4.3 for zero initial conditions, t is the natural 
period of the vehicle system at a given load. ihe fluid 
heights considered are corresponding to 1, 1/2 and 1/4 
height of the container, ihe natural frequency of the 



tn = 0.624 sec 



FIG . 4.-3 TIME Vs. MEAN SQUARE RESPONSE OF THE 
C.G. OF THE VEHICLE MOVING WITH CONSTANT 
ACCELERATION (0-868 5555 m/sec 2 ) 




t 0 = Vehicle natural period 
a = Container length 
h = Fluid height = 2.675 meter (for full) 
Vehicle acceleration = 0.86805555 m/sec 2 
t 0 = 0.624 sec { full ) 

= 0.464 sec (~ full) 

= 0.357 sec (4-full) 


FIG. 4. 4 TIME Vs. R.M.S. RESPONSE OF FREE SURFACE 
FOR DISPLACEMENT. 
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vehicle, when the fluid mass is assumed solidified, is given 
in fable 4.2. As expected., the vehicle response is a narrow 
band random process. Prom Pig. 4.3, it appears that the 
characteristics of the vehicle response do not change 
appreciably with fluid height, this is indicative of 
negligible coupling between vehicle and pendulum response. 

4.5-5 Pree surface response 

fhe r.m.s. response of the pendulum motion is the 
representative of the free surface elevation of the fluid in 
the container and is presented in Pig. 4.4* The abscissa is 
the non-dimensional parameter of time, t/t q1 where t is 
same as discussed, before, ihe ordinate measures the relative 
magnitude of the liquid surface wave. It is represented 
through the nondimens ional parameter, radian (a/2h), where, 
a, is the length of the container and, h, is the fluid height. 
The response appears to be sinusoidal. It is seen that the 
maximum relative wave height increases 1th decrease in 
fluid height,. h, in the container. As the damping is small, 
the maximum response does not decay appreciably over time. 

The steady state response occurs quite early. Numerically 
the maximum response is small enough to approximate the 
pendulum as linear. Although this aspeco is not considered, 
here but will be used in later analysis 
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Table 4.1 

Data for 2 Axl e Vehicle CR-28251 
Drawing No. SK-71605, Ref. [42] 

1. Length 7.89 meters 

2. Breadth 2.946 ” 

3- Height 2.675 " 

4- Wheel-Base 4.9 n 

5. Total sprang designed load =36.84 tons 

6. Total sprung weighr with no load =7.2 tons 

7. Lamping ratio, p = 0.2 

8. Stiffness on each wheel axle 0.358 tons/mm 

9- Moment of inertia about the 2 

empty vehicle c.g. 4245.1208 Kg-mass-meter 

Only necessary data have been listed in the table 
for the present and subsequent analysis. 


Note : 
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CHAPl'ER 5 


I'he dynamic response or a vehicle carrying 
liquid cargo on an uneven surface 


5 . 1 In troduction 

In the preceding chapter the effect of liquid sloshing 
on vehicle dynamics was investigated for an idealised model 
restricted to heave motion and subjected to a single point 
excitation due to track unevenness, i'he response was deter- 
mined by numerical integration of the non-linear equations 
of motion to simulated track unevenness. It was shown that 
for such a model, effect of sloshing on vehicle dynamics is 
not significant. In this chapter, a more realistic model 
incorporating heave and pitch degrees of freedom with two 
point excitations is considered. Ihe liquid is replaced by 
an equivalent mechanical model consisting of a fixed mass 
and a pendulum representing sloshing motion, i'he equations 
of morion of the system are non-linear, ihe track uneven- 
ness is assumed to be a homogeneous random process in space 
co-ordinate along the track. For a vehicle moving on the 
track with constant velocity, the excitation is stationary. 
For a vehicle moving with variable velocity, the excitation 
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becomes non- station ary. ihe equations of motion are linea- 
rized using the stochastic linearization technique. The 
linearized equations are transformed into a system of first 
order equations in 2n— dimensional space. ‘.Through a cano- 
nical transformation, these equations are reduced to a set 
of 2n uncoupled first order equations . ‘Ihe equations are 
now transformed from time to space domain and response 
statistics, are obtained in closed form for vehicles moving 
with constant velocity. For vehicles moving with variable 
velocity the excitation and the response are non-st at ionary 
and the associated matrices are unsymmetric. These equations 
are transformed from time to space domain. The input 
process is now homogeneous, but the co-effici ent s are func- 
tions of space co-ordinate. The response statistics are 
determined numerically. The response of the vehicle is 
also determined treating the liquid to be frozen and the 
two responses are compared to determine the effect of 
Sj.oshing of liquid on vehicle dynamics. 

5 . 2 Description of the Moving Vehicle System 

A two axle railway wagon carrying a rigid rectangular 
container filled with liquid is considered as shown in Fig. 
5.1. The vehicle moves horizontally over the track and is 
subjected to two point ground excitations. The moving vehicle 
is assumed to have heave and pitch degrees of freedom about 
its c.g. 




Loaded vehicle c.g 


-Empty vehicle c.g 


FIG. 5.1 A CONTAINER MOUNTED ON A TWO-AXLE VEHICLE 
WITHOUT TYRE, MOVING OVER A SURFACE. 
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The sloshing of liquid in the container is considered 
through its equivalent mechanical model. The liquid chosen 
is water. It is assumed that the fluid oscillations are 
small under assumed container excitations. The container 
undergoes vertical, horizontal and pitching type of excita- 
tions. Thus the parameters of mechanical model are deter- 
mined through equations (3.20) to (3-32). As the fundamental 
mode is more important than others [22], only one moving 
pendulum mass is considered sufficient for the present 
analysis. Thus the mechanical model for liquid sloshing 
has a fixed mass, m Qv , with its moment of inertia, I , and 
an oscillating mass, m , attached to the pendulum rod with 
length, 1^ . The pendulum is attached to a damper having 
damping constant, c . 

The vehicle response is determined for the reference 
point which is the loaded vehicle c.g. rhe mass, m » of 
the vehicle without container is located at a distance, , 
and has moment of inertia, In, about the pitch axis passing 
through the reference point. The c.g. of the container has 
horizontal and vertical eccentricities e-, and e 2 respectively. 
The vehicle is mounted on linear springs and dampers located 
at the corresponding axles. The spring constants are k-^ 
kg and damping constants c^, c^ at the fore and aft axles 
respectively. The whole vehicle system moves in positively 
increasing direction of x over an uneven surface as depicted 
in Pig. 5.1. 
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5 • 3 Equations of Motion 

ihe general equations of motion are determined through 
Lagrange's equations, ihe loaded vehicle c.g. has vertical 
displacement, z^, and. rotation, z^. i'he pendulum has rota- 
tional displacement , z^ measured from i~cs rest position as 
shown in Fig .5.1. 

ihe total kinetic energy (k.e.), i, and potential 
energy (p.e.), U, of the vehicle system is given hy 


2 = \ m s 2 + 7= (m z 2 + I z 2 ) + - I z 2 


.2s . 1 


2 e 


2 ' e 1 r 1 2 2 Iv 2 


+ 2 m ov t ( V r ov z 2 cos(0 ov +z 2 )) + (s-r Qv z 2 sin(0 Qv +z / 


+ I' r ' v [ ( ^l +r v^2 cos ( e v +z 2^ + 1 17 L sin z j) 2 


z^ cos z^ - r v z 2 Sin (e v +z 2 ) ) ] 


(5.1) 


u - 2 k l^ Z l + 1 1 Z 2 y l^ + 2 k 2^ Z l 1 2 Z 2 y 2^ 

+ m g [r (sin (0 +z_) - sin 0 ) + 1 (l-cos z..,)] 
v & L v v 2 v v y J 

+ m g r (sin (0 +z 0 )-sin 0 ) + m g r [-sin( 6 +z 0 ) 

ov & on v : ov 2 ov y e c e L e 2 


+ sin e e ] 


(5.2) 


where s is the displacement of the loaded vehicle c.g. at 
any time t in horizontal direction. 

The equations of motion for the system can be repre- 
sented through Lagrange's equations [25] as 
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§_ f 3 l 1 \ _ yj_ 3 u _ ... . , 

dt 3z^ 3 z^ + 3z. ~ Sji ~ ^Ai’ 1 - 1,2,3- (5*3) 

where are externally applied forces anc 1 Q Di are damping 

forces. Using (5.3), (5.1) and (5.2) one writes the equa- 
tions of motion as follows. 


(m e +m ov +m v )z l + !> OV r ov cos(0 ov +z 2 ) + Vv C0S ^ 0 v +z 2 ) ^2 
+ m v l v sin z 3 -z 3 + m v l v cos z^z? - [m^r^ cos ( ^ QY +^ 2 ) 

+ m v r v cos (\+z 0 )]'zl + (c 1+ c 0 )z n + ( c., 1, -c„lj z. 


v‘“2 /J “2 ' "'l‘~2'“l ' "'1*1 ~ 2 * 2'~2 


+ (k 1 +k 2 ) z 1 + (kjlj-kgipzj = k 1 y- L +k 2 y 2 +o 1 y 1 +c 2 j' 2 (5-4) 


[ m ov r ov cos {6 ov +z 2 ) + Vv cos ( e v +z 2 ) ]' i l 

tlI l +I o T *ov r oV 13 /( ) ' z 2 + ^v 1 v r v sin(e v +z 2- z 3» z 3 

+(m vVv 008 (8 v +z 2- z 3 ))i 3 -)(" 0T r 0 / ln (e ov +z 2 ) f 

| ■+ %i &$AA (Gv+^a^ j 

+( C 1 l 1 -c 2 l 2 )z 1 + (c 1 l^+c 2 lp z 2 + (k 1 l 1 -k 2 l 2 )z 1 
+(k 1 l^+k 2 lpz 2 + g[m v r v cos(e v +z 2 ) + m^r^ cos ( e ov +z 2 ) 


-m r cos(0 o +z 9 )] = k 1 l 1 y 1 - k ? l ? y ? + c 1 l 1 y 1 - c p l ? y 


2 " 


e 2 


1 1^1 " X 2 X 2 J 2 _r "1 X 1 J 1 ~ 2 21 


(5.5) 


Vv z 3 + VVv sin z 3 + KVv sin ( V z 2" z 3 ))z 2 


~(m v l v r v cos ( 9 v +z 2 -z 3 ))z 2 + c^ + nl^g sin z 3 = 0 

+ -WyJv £ CBS z 3 (5.6) 

If the c.g. of the mechanical model is assumed to be the 
same as that of the liquid in the container and z 2 as small 
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then one can approximate 

g[m T r v oos(e v +z 2 ) + m ov r ov c°s(9 017 +z 2 ) - m e r e eos-(B e +z 2 )] 

“ ' e 'A V 

Equations (5-4) to (5.6) are non-linear and coupled, 'these 
equations have nonlinearities in displa, cement , velocity and 
also in inertial terms. An analytical solution of these 
equations is not feasible. Hence the following stochastic 
linearization technique is used to linearize these equations. 

5 • 4 Stochastic Linearization of Equations of Motion 

There are deterministic and statistical methods to 
linearize the equations of motion. Each has its own advan- 
tages and disadvantages. These are discussed in detail in 
Ref. [4]. for the random excitation, statistical or 
stochastic linearization is preferred. Broadly speaking 
two methods of stochastic linearization are available. The 
first one was suggested by Iwan and 'fang [ 26 ] and the other 
by Atalik and Utku [4]. In the former the nonlinearities 
in damping and stiffness are considered but inertial non- 
linearities are excluded. In the latter [4], nonlinearities 
of inertial terms are also included and hence it is more 
powerful. Accordingly, in the present analysis, the method 
of Ref. [4j is adapted. 
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5.4.1 l 1 he stochasxic lineariz at ion met ho d (Ref . 


ihe nonlinear multi-degree-of -freedom system 
g (2, 'I, Z) = f(t) 


(5.7) 


is linearized to 


[m] 2 + [c] 2 + [k] 2 


?(t) 


(5.8) 


where, 

Z is the generalized displacement vector, 

g^(2,Z,2) is the total internal force acting in the ith 


?(t) 


degree-of-freedom direction., 

is the stationary Gaussian random excitation 
vector, with zero mean-, 


[mj,[c],[k] are mass, damping, stiffness matrices, ihe 
elements of these matrices are given by 

3 g'i- 

nj = E 


c = i [^] 
ag i 

= E 


i,j — 1,2,3; 


(5.9) 


where £i[ * ] is the expectation operator. 


5-4.2 


.i cation to the -present analysis 


Nonlinear equations (5-4) to (5.6) are linearized 


using relations (5*7) to (5*9). It is assumed that z 2 and 
z, are small quantities so that sin z 2 ~ z^-, sin z^ - z^. 
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cos £ 2 ~ l*0 y cos z^-1.0* The squares and products of 22 
and z are neglected* It is assumed that the excitation 
process 9 Y(s) is Gaussian in nature wit:: zero mean* Hence 
z ± are Gaussian in nature with 

B [ z i ] = 0, i = 1,2,3 (5.10) 


Using the properties of Gaussian process [32], the assump- 
tions cited above and (5.10), one gets the linearized 
equations as 

[ m i j - 1 { V + [c ip { t } + £ k id {2 i> = Uj} , i,j=l,2,3 


(5.11) 


where 


“ll = (m e + m ov + m P’ m 12 = < m ov +m v )e l’ E 13=° 


m 2 i=m 12 , m 22 = + (e^e^) 

2 


m 31 ~ m 13 J m 32 “ m 23’ m 33 “ Vh 


11 ~ ( c 1 +c 2 )» c 12 “ ^ C 1 1 1~ C 2 1 2^ 5 C 1'3 = 0 

= 0 


C 21 ~ C 12 5 C 22 


2 2 
C l 1 i + C 2 X 2 » 


'23 


c jj- — 0 , 


0, C £ 

32 33 v 


k ll = k i +k 2 y k 12 = k l' L l~ k 2 1 2 > 


1C-, 'X r\-\ i — I 


13 


'21 “ 12 s 


k 22 = Vi + k 2 X 2 ~ gffi A " (m ov +m v ) e l E [ V 2 ^ 

- Vv (e 2 +e 4 ) B f z 2 ’fcoVK® ’ 


•23 


m 


'32 


k 


33 


n (e 2 +e 4^ E t z 2 Z 3^ y k 51 = 0 
Vv < e 2 +e 4 } S [z 2 Z 3 ]s 
Vv g - Vv ^ e 2 +e 4 } 1 [ z 3 z 2 h 
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f l k l y l + k 2 y 2 + C l y l + c 2 y 2 

f 2 = k l 1 l y l " k 2 X 2 y 2 + C l X l y l - c 2 1 2 y 2 
f ^ = 0 . 

For stationary processes z, ,z it can be proved that [37] 

— J 

B [z ± zF\ = E [z i ‘z d ] (5.12) 

thus for stationary process z^ , the mass, damping and stiff- 
ness matrices are symmetric. 


5 • 5 Uncoupling the Equations of Motion for the Vehicle 
Moving with Constant Telocity 

A general method to uncouple damped linear system is 
reported in the literature [20,25,34,51]. This method, 
suggested by Boss, is specially useful for systems with non- 
proportional damping and is used here. 


Let the generalized co-ordinates are x , n=l,2,...6 


with 


6 

L 

n=l 


z . = Y u. x, z. = y a u . x , 1 = 1,2,3; 

r L . m n ' i L n m n 1 

n=l 


(5.13) 


where a n are the eigenvalues and u^ n are the elements of 
the modal column matrix 


[««] 

• ♦ * i 

[u] 

of the 6x6 dynamic matrix 

! [oj 

[ D ] = | ... ,_1 


[U] = 


(5.14) 


[I] 


-[k] _1 [m] -[k] 1 [c] ! 


(5.15) 
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where [k], [m], [cj are known from (5.11). ihe eigenvalues 
a n are < ^ e '* :erm; ^ ne( i such that 


j [D] - X [I] |=0, and X = ’ i 

Eigenvalues and eigenvectors occur in 
For a stable system, eigenvalues have 


( 5 . 16 ) 

complex conjugate pair, 
negative real part. 


and 


If , {<j> 2 • •• (bg} are the eigenvectors then 

[u] = [ (4 >1 >{4» 2 > ... {<f> 6 }] (5.17) 

[au] = [ a 1 {c{) 1 } a 2 ( cf> 2 > ... a g { 4> ] (5.18) 


are 3x6 rectangular matrices. It is assumed that [mj, [cj, 
[k] are symmetric and hence it can be proved that the modes 
are orthogonal [20,25,34]. '.Thus in generalized co-ordinates, 
the decoupled equations are [20,25,34] 

- c^x^ = 1 = 1,2,... 6 (5.19) 


where , 


B.. 


= 2a ± { cf> ± > [m] {<J> ± } +U i >' J} [c] UA (5.20) 


F. 

x 


= {4i> ( f(t)> 

= ( u 1± + l-i u 0 _. ) ( k-j y-^ +c-j y -j ) + ( u : -j — ) ( k ? y p”^Cpy p) 


1 2i M rl 11 


i 2 2i 2 2 2 J 2' 

(5.21) 


5.5.1 Equations in space-domain 

I he decoupled equations of motion (5.19) in space 


domain become 
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ff v * * * * x i-‘ x i x i = it* + 


i> + c i if Y 'i s+ 1 f ) 


+ ^ U li - 1 2 U 2i ) ( k 2 Y < S - 1 2 ) + c 2 ff T’( 8 -l 2 ))], 


where x^ 


Y » _ dY 
ds 


( 5 . 22 ) 


5 • 6 Vehicle Motion 

If s(t) is the distance covered by the vehicle along the 
track at time, t, then for constant velocity run one has 


s + v t 
0 0 


(5.23) 


where s Q is the initial distance at t=0 and v Q is the 
constant velocity. Thus 


ds • 

rrr = S = v 

dt o 


(5.24) 


5.7 Solution of Equations of Motion 


Eor constant velocity run, using (5-24), (5 >22) reduces 


v x.'-a.x. = [ ( u, . +1, u_. ) (k-, Y( s+1, ) + c n v Y 1 ( s+l-, ) ) 

oi ii R. LV li 1 2i 1 1 io r’ 


+ (u li -l 2 u 2i )(k 2 Y(s-l 9 ) + c 9 v n Y’(s-l ? ))] (5.25) 


General solution of (5.25) is the sue of homogeneous and non- 


homogeneous parts. The homogeneous part is 


B ± s 


a . s/v 
1 ' o 


( 5 . 26 ^ 


where B ± are constants and are determined from initial condi- 
tions . The nonhomo gene o us part in spatial frequency (q ) 


domain is given by 
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then the constants of integration are given by 

r 


{B} = [U] 


-1 


"1 


v, ; 


; 


(5.36) 


l 3 


where p represents mean value. 

5.7.1 Second order statistics of response 

i'he mean values for displacement, first and second 
derivatives are evaluated, using (5.31) to (5.34), as follows 

6 a s/v 6 oo .ft 

= E 0 ± U)] = l u ln B n e n + I u in i. V s1 > e ' 

i n=l n=l -oo 


B[dP(n)] = u in B n e 


a s/v 
n ' o 


n=l 


Similarly 
1 


'z. 

1 

6 

" nix 

a n 

u. B 
m n 

a s/v 
n 1 o 
e 


6 

’ nix 

2 

a 

a. 

a s /v 

-D ^ 11 ' 

e 

z . 

1 

n 

in 

n 


i = 1,2,3 


(5.37) 


(5.38) 


(5.39) 


I'he response, and its derivatives, given by Eqs.(5*32) 
to (5.34), can be termed as oscillatory processes [39, 
section 3]. Hence, after transforming to zero mean, the 
evolutionary spectral density and the covariance function 
can be determined as suggested in Ref. [39, section 4] • Ihus 
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9 


Z i Z £ (fi) =* YI ^) tN iD H^)][I_u to H a ( fl )f 


P 


=1 X P P 


n=i 


*YY^ a ) ^ H (fi) H*(0) 


(5.40) 


p=l n=l X P * n 'P 

where (*) denotes complex conjugate. Since the eigenvalues 
and eigenvectors occur in complex conjugate pair hence (5.40) 
can be written in a simpler form as 


♦z.z/ n > 
x £ 

6 

6 


= f rT fc) l 

P = 

T u. u H (a) H (-Q) 

1 n=l X P -®i P n 

(5.41) 

Similarly, 




$• • (n) 

z.z 

X £ 

6 

6 


H 
t- — i i! 

£ 

e- 

II 

V a a u. u H (q ) h (-n) 
n lp p n ip £n p'“ y Tx v ' 

(5.42) 


6 

6 9 


*z.z (!!) 

X £ 

= $ Tf (a) I 

P=1 

J^pVipVi H p (n > V-«) 

(5.45) 



i ? = 1 9 2 9 3 * 



The covariance function, for s-^ s^? is similarly expressed as 

oa 6 6 

I 

p=l 


z i z £. 


jn(s 1 -s 0 ) 


(s 1? s 2 ) = j $ YY (fl)[ l U H (Q)]r l u^jjlfe""' 1 2 

-co p=l * * -1=1 

6 6 co 

= L U ip^ n J "p ^ ^ H n (-^) 

jO — JL H—J- co 

jQ( s 1 — s _) 

e - 1 - 2 dn (5.44) 

The complex conjugate operator has been dropped on the 
argument a s given for (5.41). Similar expressions for 
covariance functions of velocity and acceleration can be 


obtained. If 



9b 


o° 3 q ( s s ) 

^pn^ s l ,S 2^ = ^ H (0) H ( -ft ) e. ^-2 ^ ? (5-45) 

— oo 4r 


then one has for s^> 

6 6 

< z.z„ = Y y u n . I 

p=l n=l lp ^n P- 0, 


(5-46) 


6 6 

K M. = Y y a a u. u. I 
1 £ p=l nil p 11 ip P n 


(5.47) 


6 6 


glpn^ s l jS 2^ 


K-i.z = v 2 y V a a u u 

1 * 0 pil nil p n ip Z* 3 s x 3 s 2 


(5.48) 


K z,z 
1 £ 


6 6 


pii ii u « 


(5.49) 


Evaluatio n of I pn ( s p> S 2 ^ and its derivatives 

On substituting for $ y ^(^) f rom (4.2) and for H (ft) 
from (5.27) into (5.45) one gets 


V = ^ iJWp +1 i u 2p )( V^ c iV e 


+ '“lp'hV 


(k 2 + jftc 2 v Q )e ' " *] [(u ln +l 1 u 2 n )(k 1 -jftc 1 v 0 )e 
3 ftlp. 3 ^( s i“ s 2 } 

(k ? - 3 ftc v ) e J — p — p 7 ■ — 7T 

" (y 2 +ft 2 )( 3^ 0 - a p )(^^o+ Q C'a) 


-jftl- 


"2 J 2 o 1 


u ln - 1 2 U 2 n^ 


(5.50) 


IE he improper integral is evaluated by the contour integration 
[30] using residue theorem. As the eigenvalue has negative 
real part, there are two poles in the upper half plane. 

One is at ft= jy and the other at ft = «p/J v 0 - Both are 
simple poles. 1 'hus finally one gets for s^ s 2 
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= f ji . n \ 

‘pn R R l %n 


p n 


pn pn' 


(5.51) 


where. 


? = ( A p Ii l +B p I) 2 ) ( VVVV ~ y(s 1~ S 2 ) 

'pn 2y "Ta p +v Q y ) ( a n -v Q y ) e 


S 1 - S 2 


G__ - + Vp> v 

J o 


pn 


u p - v o 7 2 )(“ p « n ) 


o 


D 1 = 

-yl-, 

(t rV T o) e 

yl 

D 2 = ^ k 2" c 2 Yv o^ e 2 ’ 

D 3 = 

Y 1 ] 

(k 1 +c 1 yv 0 ) e , 

P 4 = (k 2 +c 2 yv 0 )e Y 2 

"p = 

(u lp +1 lV ■ 

B P = (u lp" 1 2 U 2p ) ’ 

= 

(u ln +1 l u 2n ) - 

B n = < u ln- 1 2 U 2P ) ’ 


a In /v„ 
(k 1 +a p c 1 )e p , 

-a 1 / 

W p = (k 2«p c 2 )e P 

II 

-a 1, /v 

<Wi )e p 

a l_/v 

• V ( V“p°2 )e p ‘ 


5*7.2 Limiting cases 
5 . 7 . 2 . 1 Small vehicle velocity 


n ■ T 2 tt;C 

v7 o V = E p E n ;3 ' "> 


lira (Jt„ + G„J. , 
0 


,- y(s 1 - S 2 ) 

V = 2 tlT^ 6 


o 


lim G ^ = — r 
0 P n - 2 


P n 

1 n--„ _ „ a p^ S l" S 2 ' )//v o ri , 2, 1 .2 „ 2 / „ 2 , „ 2 > 


V 


a p ( “p' t V V° 


lim v Q e 


[k£+k 2 -a^(c£+c|) 


a L/v -a L/v 

+ ( k l +a p 0 l><V“ p° 2 > e P + (ki- a p c l )(k 2 + 0 : p c 2 )e P J 


where L = 1-^ + 1_ . 


(5.52) 
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is a p has negative real part 


lira v e 
0 


a p< L+S 1“ S 2> / v 0 


“0 , for all s-j^ s . 


, . -a(l-s+s)/ v 

limve^ d °-*-n 

v o -0 0 ? 


if h<(s 1 -s 2 ) , 


" » if 1 > (s 1 -s 2 ) 5 

as the exponential term increases much faster than v c 
reduces. Hence 

lim G ->oo 
v-0 pn 


if 1 > (s 1 -s 2 ) 


lim 

vi' 0 
o 


'pn 


(5.53) 


ior mean square response s p =S 2 " i'his implies, from (5*46) to 
(5-48), that the mean square response in the neighbourhood of 
zero velocity is large. 

5 . 7 . 2 . 2 L arge vehicle velocity 

(i) lim I = lim ( F + a \ 

v pfl HA ^ pn pn- 
0 1 0 

Looking to the expressions of tf and G pn , from (5-51), one 
can easily prove that 

.2.„ t*> .2 


Wl +B P B n c 2 + < Vi 

2y 


-yL 


+ B A e Yl ) c n c. 

lim j? - p u ~ b n z ~ p n p r 1 2 

Y P^* 

0 


(5.54) 


lim G- = 0 

Y -> i00 

O 


=IT^7 [WlH B n 0 f«l c 2 <Vn S " Tl+B pV Y 5] 


(5.55) 
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i‘his implies, from (5.46) and (5.47), that at large velocity 
the mean square response for displacement and velocity tends 
to he independent of velocity. 


(ii) lim v 


2 3I S-=ltm ^ r..2.2 


v 0 -*» ° 3s l 3s 2 v 0 -vco ^n 


p p [y^vfF G ] (5.56) 

±LR„ 0 pn p prx J 


Looking to the limits evaluated at (5.54), or directly one 
finds that this limit tends to « . This implies, from (5-48) 
that the mean square response for acceleration tends to be 
large at large velocity of travel. 


5 . 8 Shear-Force and Bending Moment Response 

Shear force and bending moment coming on the bottom 
corner of right container wall are determined assuming the 
wall to be a cantilever. 1'he forces are assumed to act 
along the vertical centre line of the container wall which 
is fixed at the bottom. From Fig. 5.1, the shear force 
is given by 


where , 
F 


F 


ov 


F = F + F , 
s v ov 


| "V ft (V 3 cos - V 2 8111 (e v +z 2 }) 

1 xr (~- - r z _ sin (6 +z 0 )) 

2 ov dt v dt ov 2 ov 2 


(5-57) 


As Zg z 3 are quantities one has 

?s = b o + b l z 2 + b 2 h + VV 2 


(5.58) 
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where , 


1. 


b o = I“oy 3 ' ' b l = - 3 K r Y 8iJ1 e „ + =>„„r,„sin 9 „ J 


V 


0 ? ov 


h 2 - 2 W b 3 “ ' 2 (m v Wi ov ) e l 
Similarly the bending moment is given by 


M b = (§ + e } r 


<§ - b 5> J oy 
where, h is the height of the liquid . 


ov 


(5.59) 


(5 - 60 ) 


Or, M b 
where, 
r o 


r Q + r^z ^ + r 2 z^ + r^z 2 z 2 


(f - e 5 )b o’ 

- | [ m v (|+® 4 )V ln 6 Y + m 0Y ( f- e 5 )r 0Y Sin 0 OYi' (5 ' 62) 


( 5 . 61 ) 


r~ = 


( 2 +e 4 )b 2 , 


r 3 “ 


[mv(f^ e 4 ) + m n jj-e^)]e-, . 


ov 2 5 


5*8.1 Mean Value response 


As the mean of z. is zero hence 




= + N B[z 2 z 2 ] . 


0 


(5.63) 


for a stationary process z , one has 
E[z 2 z 2 ] = -B[z|] • 

This is a second order term and is neglected. 


( 5.64) 




= b. 


0, for constant vehicle velocity. 


(5.65) 


Similarly, 
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— 0 f for constant vehicle velocity. 


5*8*2 Mean square response 


(5*66) 


Under Gaussian assumption one gets, after neglecting 
the second order terms, 


2[(F S -^) 2 ] = b 2 E[z 2 ] + b 2 E[ z^J + gb^E [z^] .(5 .67) 
Similarly, 

E[(iV^ )2 ] = ^S[z 2 ] + r 2 E[z 2 ] + 2r 1 r 2 E[z 2 i 3 J . (5*68) 

5 • 9 Equations of Motion for the Vehicle System with 
Solidified Mass ~ ’ ‘ ~ 


The equations of motion for the vehicle system carrying 

equivalent solid cargo mass, M , having moment of inertia, 

I Q , can be written as follows, 
sy 

[m] { z^} + [c]{ z.} + [k]{ z J = (f J , i = 1,2 (5*69) 


where 


[m] : 


[k] : 


m +M 
e o 


0 


r 


c-, +c , 


L° 


I l +I sy +M b (e l +e 2 ) 


2 2 ' * 9 t : 


Cl 1 ! c 2 1 2 


1 


1 1 C 1 -1 2 C 2 C 1 1 1 +C 2 1 2 


L 


k l +k 2 1 l k l -1 2 k 2 
1 l k l~ 1 2 k 2 k i 1 i +k 2 1 2 


f 


1 


| k l y l + k 2 y 2 + C l y l + C 2 y 2 


J i 2 _ 


l k l 1 l y l- k 2 1 2 y 2 +C l 1 l y l- C 2 1 2 y 2 
These equations are linear and coupled with 2 degrees of 


freedom. The solution of these equations and response 
statistics can be obtained as pointed out in Sections (5*5) 
to (5*7). . 
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5 • 10 Vehicle Moving with Const ant Acceleration 

1'he linearized equations of motion are given by (5-11)* 
For the variable velocity case the response becomes non- 
st a/cionary. In equation (5.11), the stiffness matrix contains 
expectation terms which are nonstationary . Hence stiffness 
matrix is unsymmetric and the equations of motion cannot be 
uncoupled by the procedure outlined in Section 5*5. the 
equation (5.11) can be transformed to 

{q} + [A] (q> = [M]" 1 (P(t)} 


(5.70) 


where 


[A] = [M] _1 [K] , 

m = 


f{z}] 

'{0} ] 

l, { F(t)} = J 

. , , i 

{flt)} 

» J 

\{z}\ \ 

V J 


[0] [m] 


-[m] 

[0] 

, [m] [c] 

, [K] = 

[0] 

W 


[M] = 

L m ]» [c], [k] are ilxN square matrices whereas [14] , [K], [A] 
are 2Nx2N square matrices. {q(o)> and [ Kt)} are vectors 
with 21 elements, net 

{q} = [D] (*> (5 ’ 71) 

On substituting (5.71) into (5-70) a no! premultiplying by 
[u ] -1 one obtains 

{i} _ [a] {x} = (g(t)> 


(5.72) 
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where, (g(t)} = U" 1 [M] -1 {E(t)} (5.73) 

[a] == -[U]" 1 [ A] [ U] (5.74) 

I'he matrix [a] is a diagonal matrix if [ u] is the modal 
matrix of [A]. 1'hus the equation (5*72) represents a system 
of 2I\I first order uncoupled equations in time domain. 


5 .10.1 Solution of Uncoupled equations in space co-ordinate 
Equations (5*72) in space co-ordinate reduce to 
v x^ - omx^ = Gh ( s ) , i = 1,2,.. 2u (5.75) 


where 


v = ft = ^o +2C o s ’ = 

and v Q is the initial velocity, C 
tion. G^(s) is the corresponding 
transformed in space co-ordinate, 
of (5-75) is given by 



ds 


is the constant 
o 

element of g^(t) 
The homogeneous 


(5.76) 

accelera- 

when 

solution 


x ih = B i eXp ( a il/fvf+2C^/C 0 ) (5.77) 

where are constants of integration. 

The nonhomo geneo us solution is expressed in Duhamel 
integral form and is given by 

x. = J h.(s,p) G.(p) dp (5-78) 

m 0 i 

where the impulse response function h^s,p) is given by 


h ± (s,p) 


V ( v q 2 C 0 P) 


exp 


[ ct i ( ¥ ( 


-o + 2C o s) - fK + 2C o p))/G o ] 


(5.79) 



(5.80) 


thus "tli© complete solution is given by 
X. (s) = X + X. 

i in m 

Finally using the relations (5.70), (5.71), (5.80), one can 
obtain the response z^. 


5 • 11 Results and Discussions 

i'he vehicle considered, for the present study is the 
railway wagon whose data are given in ia'ble 4.1. i‘he fore 
and aft axle stiffness and damping ratio are same, ‘i'he 
container is filled with water and is mounted on the vehicle, 
ihe response of the vehicle is determined for liquid cargo 
and also for the equivalent solidified mass, i'he liquid is 
replaced by an equivalent mechanical model whose parameters 
are given by Sqs . (3.20) to (3-32). Only the fundamental 
mode of liquid sloshing is considered as the contribution to 
the total wall forces and vehicle dynamics from higher modes 
is expected to be small [l,22j. In th-- present analysis, the 
fixed mass location of the mechanical model is determined 
such that the c.g. of the model is the same as that of the 
liquid in the container. I'he base excitation to the vehicle 
is provided by the vertical track ro ugliness whose p.s,d. is 
given by Eq. (4.2). 

5.11.1 Vehicle moving with constant velocity 


I'he vehicle response is studied for nine different 
configurations of the container. 1‘he c.g. of the empty 



vehicle is taken as the origin ( 0 , 0 ). rhe c,g. of the 


container is referred to the origin as (x,y), where the . 


first co-ordinate refers to the horizon'll distance and the 

, a 

second co-ordinate refers to the vertical distance. For 

example, (a^/ 4 , h^/ 4 ) refers to a container with its c.g. 
at a distance of a -^/4 to the right of the origin and at a 


distance 

of h^/4 above 

the origin. 

ihe nine cases consi 

dered are 

given below: 



Case 

Length 

Liquid-depth 

Container c.g. 

1 

a l 

h 

( 0 , 0 ) 

2 

? t 

V 2 

( 0 , -1^/4) 

3 

»! 

V 4 

( 0 , -32^/8) 

4 

a x /2 

V 2 

(a 1 /4 s 0 ) 

5 

It 

V 4 

(a 1 / 4 J h-^/ 8 ) 

6 

It 

V 2 

(a 1 /4, h x /4) 

7 ■ 


V 4 

(a 1 / 4 , f^/ 8 ) 

8 

It 

V 2 

( a 1 /4, 2^/4) 

9 

It 

V 4 

( a 1 /4, 3h-j_/8) 


where a.^ = 7.89 meter, h^ = 2.675 merer. 

It has been seen in Sec. 5.7.2 that the response 
tends to a lange value in the neighbourhood of zero vehicle 
velocity due to a mathematical singularity. (Therefore, 
the velocity of interest has been chosen beyond 5 m/sec. 

Ihe upper limit of the velocity is Kept ax 95 m/sec as it 
is considered sufficiently large for the vehicle of interest. 
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It erative procedure 

It is seen from Eq. (5.11) that the elements of mass 
and damping matrices are constant. Some of the elements of 
the stiffness matrix include expectation terms. I'hese terms 
are not known initially. Hence, for response calculations, 
following iterative steps are undertaken to determine the 
expectation terms involved: 

(i) Initially the expectation terms are assumed to be zero 

(ii) the response is determined and the terms involving 
expectations are evaluated; 

(iii) with the expectation terms of step (ii) , the response 
is determined and also the expectation quantities; 

(iv) expectation quantities of steps (ii) and (iii) are 
compared to establish the convergence; 

( v) if the convergence is not achieved, within specified 
limits, steps (ii) to (iv) are repeated. 

For all the vehicle configurations, it was found suffi- 
cient to keep the limit of convergence within 1 percent of 
the quantities of interest and 2 to 3 iterations were found 
sufficient . 

i'he vehicle response for the case of solidified mass 
of liquid is determined using the equations of motion 
outlined in Section (5.9)* 

I'he two cases of cargo as liquid, and the same as soli- 
dified mass, are distinguished with the letters Is’ and »N* in 
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various figures to present the results, ihe superscript 's' 
reiers to the case of sloshing liquid and the superscript 
'N' refers to the case of solidified mass of the liquid, i.e. 
the case wi dh no-sloshing, in each figure, which accompanies 
a table, the numbers indicate the container location, size 
and liquid height . 

5 .11.1.1 • S . — h eave response for displacement and velocity 

In Fig. (5.2a, b), the r.m.s. displacement and velocity 
response for heave mode of the loaded vehicle c.g. are presented 
for nine different cases. It is seen that at small velocities, 
the response fluctuates rapidly for all the cases considered. 

For vehicle velocities more than 15 m/sec, the r.m.s. value 
of the displacement response increases slowly with velocity, 
f'he r.m.s. velocity response fluctuates up to 30 m/sec. It 
increases rapidly up to 60 m/sec and then increases gradually. 
Both the responses tend to be independent of vehicle velocity 
at ixs large value, this is in accordance with the limits 
discussed in Sec. 5.7.2. 

For centrally mounted container, (cases 1,2,3), the 
equations of liquid sloshing and vehicle vibration in heave 
mode are uncoupled. Hence sloshing has no effect on vehicle 
response in heave mode as is evident from Fig. (5.2a,b). 

For eccentrically mounted container the r.m.s. response 
for sloshing case is invariably less, as much as 45 to 50 $, 
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than no— sloshing case (compare cases 4 to g in Pig. 5.2a). 
lhe r.m. s. velocity response for sloshing case is greater 
chan no— sloshing case (cases 7,8) for some velocities. But 
the range of such velocities is narrow and towards lower 
end. Beyond vehicle velocity of 40 m/sec (144 km/hour), the 
response with sloshing is invariably less, as much as 45 
to 50^; than without sloshing case. For the same horizontal 
eccentricity and mass of the liquid, larger vertical eccen- 
tricity causes greater reduction in the response for sloshing 
liquid (compare the cases (4,6) and (7,9) in Pig. 5.2a,b)). 
1'his implies that for larger eccentricixy the absorption 
effect of the sloshing liquid on the heave response of the 
vehicle is significant for the cases considered. 

5.11.1.2 R.M.S. pitch response for displacement and velocity 

Bine vehicle-container configurations are chosen to 
present the r.m.s. displacement and velocity responses in 
Pig. (5.3a,b). It is seen that at small vehicle velocities, 
the response fluctuates rapidly, lhe fluctuation ceases 
around velocity of 15 to 20 m/sec= After this, the response 
increases and a peak is reached again in the velocity range 
of 20 to 40 m/sec. Por displacement response there is gradual 
reduction in response after the peak (cases 1,2, 3, 7, 9)* I'n 
some other cases the response increases gradually (cases 
4, 5, 6, 8). Por velocity response there is gradual reduction 
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In response after the peak which occurs in the range of 
15 to 45 m/sec. In some cases, the velocity response tends 
to become independent of vehicle velocity quite early 
(Fig. (5-3b), cases 4, 5,6,8) . 

For centrally mounted container, the equations of 
liquid sloshing and vehicle vibration in pitch mode are 
coupled. Sloshing brings down the pitch response by several 
orders of magnitude. With the decrease in liquid depth, 
the vehicle response in pitch reduces considerably (cases 
1,2,3 of Fig . 5 . 3a, b) . 

For eccentrically mounted containers, the sloshing of 
liquid reduces the response over most of the velocities. 

For some velocities the response increases. But the range 
of velocities is narrow and towards lower end (cases 4 to 9)* 
For the same mass of liquid and horizontal eccentricity, the 
sloshing of liquid reduces the response with increase in 
vertical eccentricity. Compare cases *.4,6) and (7,9). It 
is expected that with increase in eccentricity, there is 
more of sloshing. 'Thus it implies that more the sloshing of 
liquid, larger is the reduction in vehicle response. For 
a typical vehicle velocity of 40 m/sec (144 km/hour) , the 
reduction in displacement response is of the order of 6 to 
Y5% and that of the velocity response is as much as 30 to 
60 c /o (cases 4 to 9) • 



119 


5.11.1.3 R.M. S. free surface displacement response 

Free surface displacement response statistics can be 
used to predict the spilling of contained liquid under random 
excitation. If the second order statistics of the response 
is known, the maximum response can be predicted. 

The pendulum displacement represents the free surface 
displacement of the contained liquid. In fig. (5.4a), the 
free surface response is presented. The non-dimensional 
parameter (radian- a/2h) , represents the relative wave 
height at the free surface of the liquid (a = length of the 
container, h = depth of the liquid), The surface wave height 
decreases with velocity and. approaches a constant value as 
the vehicle velocity increases. This is in conformity with 
the limits discussed in Sec. 5.7.2. 

The surface wave height increases with decrease in 
liquid depth (cases 1,2,3? 5,6? 7,8). Por eccentrically 
mounted container the relative surface wave height is much 
larger for same length and depth of the liquid (cases 4,6$ 
7,9). In general, the eccentrically mounted container offers 
more surface wave displacement than the centrally mounted 
container for most of the velocity range. 

5.11.1.4 R.M.S. pendulum velocity respo nse 

The pendulum velocity represents cne tree surface 
velocity of the contained liquid. The response fluctuates 
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rapidly in the small velocity range (Pig. 5. 4b). ihe fluc- 
tuation ceases in the range of velocities 10 to 25 m/sec. 
After this, the response increases and reaches a peak and 
then starts decreasing. It is seen that rhe response 
approaches a constant value rather early (approx, velocity 
of 45 m/sec.) as compared to the displacement response 
(approx, velocity of 80 m/sec.). 

It is seen that the response decreases with decrease 
in liquid depth (cases 1,2,3; 5,6; 8,9). It is alsu seen 
that for same length and depth of liquid, che response 
increases with increase in eccentricity of the container 
( cases 4,6; 7,9) . 

5.11.1.5 R.M.S. shear force and bending moment response 

I'he r.m.s. shear force and. bending moment response are 
presented in Pigs. 5-5 and 5-6 respectively, i'hese are 
determined using the relations (5-67) and (5-68) for all the 
vehicle-container configurations. l‘he ootal force exerted 
by the liquid is assumed to act on tie vertical centre-line 
of the side wall, i'he shear force and bending moment are 
determined on the bottom edge of the side walls. 

It is seen that the response is larger in the small 
velocities range and. fluctuates rapidly, ine fluctuations 
may continue up to a vehicle velocity of 25 m/sec. After 
this, the response increases to reach a peak and then starts 
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decreasing. It is expected that the response will start 
increasing again (case 1). It is to he noted that at large 
velocity, the response tends to a large value in accordance 
with the limits discussed in Sec. 5.7.2.2(ii). 

Table 5.1 gives the mechanical model parameters, 
hydrostatic and r.m. s. shear force and bending moment values 
for various depths of the liquid in the container. The 
hydrostatic values may be treated as the mean values. The 
coefficient of variation of dynamic forces can be defined 
as the ratio of the r.m. s. value to its mean value. Using 
the results of Figs. 5.5 and 5.6 the co-efficient of varia- 
tion can be determined. In Table 5.1(iii), the co-efficient 
of variation for dynamic shear force and. dynamic bending 
moment are given. A typical value corresponding to the peak- 
value observed in each of the nine cases is given in Table 
5.1(iii). It is seen that the co-efficient of variation 
increases with decrease in liquid depth (cases 1,2,3? 4,5; 
5,6; 7,8; 8,9). Similarly it increases vat h the increase in 
eccentricity (cases 4,6; 7,9). Although the co-efficient 
of variation increases with decrease, in liquid depth, it 
is seen from Figs. 5-5 and 5«6 that tno absolute value of 
the dynamic force increases with increase in liquid depth 
(cases 1,2,3; 5,6; 7,8). From Table 5.l(iii), it is seen 
that the co -efficient of variation for the shear force can 
be as large as 1.556 and for bending moment it can be 0.218. 
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i'he lowest value of the co-efficient of variation, corres- 
ponding to the centrally mounted and completely full container, 
is 0.113 for shear force and 0.004 for bending moment. If 
one assumes the maximum dynamic response as three times the 
standard deviation, the maximum dynamic bending moment 
( c a.s e 1, Table 5.l(iii)) is nearly 1.2% of the static bending 
moment . Similarly the maximum dynamic shear force will be 
nearly 34$ of tne static value for the same case. Thus the 
results indicate that the effect of liquid sloshing is not 
significant in the design of the container walls for strength. 
The dynamic effect of liquid sloshing may, however, be impor- 
tant from the point of view of fatigue. 

5.11.2 V ehicle moving with constant acceleration 

Vehicles generally accelerate and decelerate during 
starting and stopping respectively. lor the dynamic analysis 
of the vehicle system under variable velocity run, equations 
of motion (5-4) to (5-6) are considered The parameters of 
vehicle and mechanical model considered aro the same as for 
constant velocity case. The vehicle maintains constant 
acceleration and is assumed to acnieve a velocity of 150 
km/hour over a distance of 1 km from rest. The distance and 
velocity at any time, t, can be determined from Els. (4.8) 
and (4.9). The vehicle rolls over a railway track whose 
spectrum is given by Eq. (4-2). The track roughness provides 
vertical excitation to the vehicle system. Two sets of 
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eq.ua.ti.ons of motion are considered for the study, i'he first 
set of equations (5.4) to (5.6) includes sloshing of liquid 
and the second set of equations (5.69) belongs to no-sloshing 
case . 

5.11.2.1 Numerical solution of the equations of motion 
Simulation 

the two point inpux excitations, y^t) and y (t), are 
the road roughness at a distance, L, apart (i = base length 
of the vehicle). Using the power spectral density of rail- 
road roughness, given by Eq. (4.2), the input excitations, 
y^(t) and ygCt), can be generated as discussed in Chapter 4- 

Numerical integration and averaging 

Equations of motion (5.4) to ( 5 - 6 ) and (5.59) are 
solved by Runge-Kutta algorithm for zero initial conditions. 

As the input excitation is random in nature ensemble average 
of sample functions are taken. 5 sample functions are 
considered for averaging. Eurther averaging is carried out' 
over five consecutive readings of the sample response function. 

Presentation of results 

Three different vehicle-container configurations are 
chosen for the constant acceleration of the vehicle. As in 
constant velocity case, the results are presented for the 
sloshing and no-sloshing cases. These are distinguished with 
superscript 's’ and 'N' respectively for sloshing and no-sloshing 
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cases. The three different configurations are numbered and 
tabulated in the various figures as in constant velocity 
cas e . 

5.11.2.2 Displacement response 
R.M.S. heave response 

It is seen from Pig. 5.7, that the response with 
sloshing and without sloshing are almost identical. This 
implies that for centrally mounted container, the heave 
response remains virtually unaffected by the sloshing of the 
liquid. The response is a narrow band random process with 
the period same as the heave period of the vehicle (Table 
5.l(ii)) . The response is also compared with the vehicle 
having single base— point excitation described in Chapter 4- 
It is seen that the vehicle response with single-point base 
excitation is in general larger than the present case with 
two-point base excitations. This is in agreement with the 
results obtained in reference [53]. 

R.M.S. pitch response 

Por the same vehicle-container configurations, a.s for 
heave case, the results are presented in Pig. 5*8* The 
results are also compared with no— sloshing case. It is seen 
that, except for a very small period of time considered, 
the response without sloshing case is larger than with- 
sloshing case. The liquid depth has a considerable effect 
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on the response statistics. Small liquid depth reduces the 
response by several orders. 1‘his is in accordance with the 
observations already made as in the constant velocity case. 

R.M.S. free surface response 

The results for this case are presented in Pig. 5 • 9 • 

The pendulum displacement is the representative of the free 
surface displacement. The non-dimensional parameter 
(radxa/2h) represents the relative liquid free surface 
displacement. It is seen that the liquid surface displacement 
reduces with the increase in liquid depth. The response is 
a narrow band random process. It has a same period as that 
of the pendulum (see Table 5.1(ii)) for the completely full 
container. Por half -full and one-fourth full containers, 
the period is same as that of the pendulum but high frequency 
of excitation is also superimposed on the response. 

One can compare the present case response, having two- 
point base excitations for the vehicle, with the correspond- 
ing free surface response for the vehicle having single 
point base excitation given in Pig. 4 «4* It is seen that 
the single source excitation system of Chapter 4, keeps the 
free surface response unaffected by the randomness of the 
"fox* roughness . It indics-t es that heave venicle model has 
weak coupling with the sloshing of .liquid. The present two- 
point base excitation model has significant coupling with , 
the sloshing of liquid. ■’ 
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5 • 11 .2.3 Acceleration response 
R.h.S. heave response 

In. Pig. 5.10 the acceleration response, in non-dimen- 
sional form, is presented. The ordinate is the ratio of 
vertical acceleration to acceleration due to gravity. As 
expected, the vertical acceleration increases with decrease 
in liquid depth. The acceleration response is a narrow 
band random process and at the start has a impulsive nature. 

The response is same for sloshing and no -sloshing cases. It 
implies that coupling effect on heave response due to 
sloshing is negligible. Further the acceleration response 
is small (0.01 to 0.07 times the acceleration due to gravity) 
for the kind of vehicle under consideration. 

-R.M.S. pitch response 

The acceleration response in pitch is shown in Fig. 

5.11. As expected, the response reduces with increase in 
liquid depth and is a narrow band random process. For most 
of the time range covered, it is larger for no-sloshing case 
than the sloshing case. It is also seen that the difference 
between these cases increases with reduction in liquid depth. 
This is similar to the one observed for constant velocity 
case. Thus sloshing of liquid has an absorption effect. 

-R.M.S. pendulum response 

The acceleration response char act eristics of the pendulum 
represents the mean free surface acceleration of the contained 
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liquid . It is seen that the acceleration of the pendulum 
increases with decrease in liquid depth. The relative 
increase in acceleration with decrease in liquid depth is 
almost of the same order. 

5 .12 Absorption Effect of the SI o shin- Liquid 

In Sec. (5.11.1) and (5.11.2) the response for constant 
velocity and acceleration of the vehicle are presented. Ihe 
response in heave and pitch mode of the vehicle are compared 
with the two cases of sloshing liquid and equivalent solidi- 
fied mass of the liquid. 

It is observed that the sloshing of liquid has the 
tendency to reduce the response of vehicle in heave and 
pitch mode of the vehicle vis-a-vis the vehicle carrying 
equivalent solidified mass of the liquid. This implies that 
the sloshing of liquid acts as an absorbing mechanism for 
the vehicle. For the given input, the energy is shared by 
the vehicle mass in motion and. the sloshing liquid mass. 

The energy to the sloshing liquid mass is transferred through 
the vehicle, thus reducing the vehicle response. In contrast 
to this phenomenon, the energy shared bj The solidified mass 
of the liquid becomes the integral part of the energy shared 
by the vehicle system. Therefore the vehicle response reduces 
due to the sloshing of liquid. 

The vehicle response reduces with the increase in slosh- 
ing of liquid. This is usually accompanied by reducing the 
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liquid depth. This observation is verified both for constant 
velocity case and constant acceleration case. It is known 


that the vibration absorption becomes more effective with a 
reduced ratio of the natural frequency to the exciting fre- 
quency [52]. A container mounted on a vehicle is subjected 
to narrow band vehicle excitation, as ciscussed in Sections 
5.11.2.3 and also in 4-5-4, with the central frequency 
corresponding to the natural frequenc 3 r of the vehicle. In 


table 5 .l( ii) , the natural periods of vibration, for the 


vehicle with solidified mass of the liq-id in centrally 
mounted container, are given for heave and pitch modes. It 
is seen from table 5-1 (see also Table 4.1) that the ratio 
of natural frequency of liquid sloshing to that of vehicle 
natural frequency reduces with decrease m liquid deptn. 
This ratio, for liquid depths^ (=2.675 meter) to 1^/4* 
varies from 0.171 to 0.0584 for pitch mode and 0.174 to 
0.0573 for heave mode of the vehicle. J - llus sloshing 
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CHAPTER 6 


SUMMARY OP RESULTS AND CONCLUSIONS 


In the preceding chapters the results of the investiga- 
tions of the following problems in the sloshing of liquid in 
moving containers have been presented*. 

(i) The stability of the plane free surface of liquid in 
a rectangular container with flexible bottom, under 
periodic vertical excitation? 

(ii) development of a mechanical model for liquid sloshing 
in a rectangular container under The simultaneous 
action of vertical;, excitation and horizontal motion, 
and 

( iii) study of the dynamics of a vehicle carrying liquid as 
cargo and comparison of the results with a vehicle 
carrying an equivalent solidified mass of the liquid 
under identical random track induced excitation. 

The results have been discussed in derail in ocxch 
chapter. In the sequel some of the important results, 
conclusions and recommendations for further investigation 


are presented. 
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(1) The stability of the plane free surface of liquid in 
a rectangular container with flexible bottom under vertical 
harmonic excitation has been investigated. Stability bounda- 
ries have been established and compared with that of container 
with rigid bottom. Since a container mounted on a vehicle is 
subjected to narrow band vehicle excitation with the central 
frequency corresponding to heave degree of freedom, tile 
results of the stability analysis can be used as a guide for 
the behavior of liquid in vehicle mounted containers. The 
work presented in the thesis may be extended to include the 
flexibility of the side walls and consider narrow band random 
excitation. 

(2) The sloshing of incompressible and inviscid liquid in 
a rigid rectangular container, under simultaneous action of 
constant horizontal acceleration and vertical periodic 
acceleration, can be represented through its equivalent 
mechanical model. The model consists of a fixed mass and 
infinite number of oscillating pendulum masses. The parameters 
of the mechanical model are independent of the excitation 
characteristics. The fixed mass location of the mechanical 
model can be decided on the assumption that the mechanical 
model c.g. is the same as that of the liquid. This is 
assumed in the present analysis. The model is valid for 
small displacement and stable free surface conditions of the 
containing liquid. It is capable of incorporating important 



143 


modal characteristics of the oscillating free surface of the 
containing liquid and, forces and moments exerted by the 
sloshing liquid on container walls. 

(3) The vehicle unevenness of the track is treated as a 
homogenous random process. For an analytical treatment of 
vehicle response, it is convenient to represent the power 
spectral density of the vertical unevenness of the track by 
the following expression; 

c 

$ Yy .(f) = , one sided p.s.d., ft 2 /( cycle/ft) ; 

11 y 2 +f 2 

where C and y are constants. It is shorn that the above 
expression fits the measured p.s.d. closely. 

The equations of motion of the vehicle-container system, 
with liquid replaced by an equivalent mechanical elements are 
non-linear, ihe vehicle response is studied, in two ways 
under constant acceleration motion. In the first the vehicle 
is assumed to respond under heave mode and, in' the second it 
is assumed to respond under heave and pitch modes. In the 
heave mode alone, the vehicle is exposed to single point base 
excitation and. in the second it is exposed to two— point base 
excitations. 1'he response of the equations of motion is 
determined numerically* using the simulated track unevenness. 

In case where acceleration is not constant, it can be 
approximated by segments of constant acceleration to deter- 
mine the response. 
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It is found that the effect of liquid sloshing on the 
vehicle response is negligible for the heave model, ihe 
response characteristics of the vehicle model with heave and 
pitch degrees of freedom are, however, significantly affected 
by the liquid sloshing, i'he sloshing of liquid reduces the 
vehicle response due to absorber effect, ihe free surface 
response in heave model remains unaffected by the randomness 
of the track unevenness whereas in the heave-pitch model, 
the free surface response has random characteristics. 

(4) I'he vehicle response with heave model is found to be 
larger than the heave-pitch model. Ihe heave model has 
single point base excitation and the heave-pitch model has 
two-point base excitations. 

( 5 ) the coupled vehicle-liquid sloshing behaviour is 
investigated for a vehicle model with heave and pitch degrees 
of freedom under constant velocity and acceleration over a 
railway track, hine different configurations of the container 
location and liquid depth are considered, ihe liquid is 
replaced by an equivalent mecnanical model. Ihe governing 
equations of motion of vehicle— container system contain 
non-linearities in inertia, damping and stiffness terms. 

i'he equations are linearized using the stochastic lineariza- 
tion technique. 'I'he response statistics are determined 
after transforming the equations to space co-ordinate. With 
space co-ordinate as independent variable, the input process 
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is homogeneous and a unified treatment is possible for deter- 
mining the response statistics both for a vehicle moving 
with constant velocity and constant acceleration. 

Some of the important results obtained from the investi- 
gation of the response of a railway wagon carrying liQ. uaa aS 
cargo are: 

(i) Sloshing has a significant effect on vehicle dynamics 
for eccentrically located containers? 

(ii) the liquid acts as a vibration absorber at vel 00 ^"^ 63 
more than 40 m/sec. ihe vehicle response is reduced 
due to transfer of energy from vehicle to liquid aS 
in the case of two degrees of freedom system. ' x ‘he 
absorber effect is more pronounced in pitch-deg ree o:f 
freedom. In some cases, for velocities less than 

40 m/sec, liquid sloshing increases the response 

(iii) the relative wave height at the free surface increases 
with decrease in liquid depth and with increase in 
eccentricity .of the container c.g. w.r.t. vehi°d e c, §' 
The liquid sloshing is found -do reduce the vehi°de 
response under certain conditions. Eonce it may be effectively 
used as a means of vibration absorption to improve the riding 
quality of vehicles. 

(6) It is found that the effect of liquid sloshing xs no " t 
significant in the design of the container walls for strength. 
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The dynamic effect of liquid sloshing may, hoover, be impor- 
taut from the point of view of fatigue. 

(7 ) For the design of the wagon under-carriage, the effect 

of liquid sloshing on vehicle dynamics may oe neglected, 
the specific case considered, the model with no-sloshing 
gives conservative results at hign velocities. 

(8) The model of the railway wagon considered m the s y 

includes only heave and pitch degrees of ireodom. 
wagon has all the six degrees of freedom and is subjected o 
track induced excitation at each wheel due to vertical and 
lateral irregularities . A more realistic analysis must t e 
in to account other degrees of freedom of the wagon, particular 

lateral and rolling, and lateral track irregularities also 

+ . + i0 n will require a more general equiv en 
Such an investigation win 

n n . for iiquid sloshing, ike investigation 

mechanical model for liquiu 

, . . +n cover the effect of various parameters, 

should he extended to cover 

, , n f vehicle suspension* 

such as the stiffness and aampm* of vchi 

on the vehicle dynamics. 
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jiPPEKDIX A 


SO MS 


& EHQEBBTIBS 0? MATHISU-HIIi SQailOS 


Following Eiathieu-Hill equation 

the stability characteristics. 


a!? +'-n 2 [i - f(t)]y = 0> 

at 2 


is considered to study 


(A-2.1) 


where 


f(t) is a periodic function with a period 


■i = 


2ri 

03 


/ h 0 i \ does not change its 
- f(t), equation U-2-U aoes 

Since f(l+t;- I v^> that if 

... . „fth e period T to t. This umpires that 
form on addition of P ^ 

.. n vft+i’) is also the solution of 
y (t) is the solution then y( 

(A-2.1). 

- r o' y\ f\ is linear in 
■ ( L ? 1) is Of s econo or an - 

SQuaxron + x , . 

H nee if y (t) and y 2 (t) are two independent so u ion 

y ‘ “' 1 . a y (t+l‘i are also its solutions, 

of (A-2.1) then y 1 (‘ t+i ) and y 2 l 

fp expressed in the linear 
xhese last two solutions can he excess 

/ + \ o-nd y (t) ns follows. 

combinations of y-^( 2 

(A-2.2-1) 


y 1 (t+T) = a u + a 12 y 2 (X) ’ 

y 2 (t+l) = a 21 y x (t) + a 22 * 2 {t) ’ 
where a, , are constants. 


U-2.2.2) 
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i 


Thus, the addition of the period S to t results in a 
..inear transformation of the initial system of solutions. 

It is to be noted that other combinatioi s of independent 
solutions are possible on the right hand side of the last 
two equations. In particular, one can try to choose solutions 

y*(t) and y*(t) such that the co-einciente a 12 and a 2 l 

1 . s . , ,n>, 0 transformation in this case 

equations U- 2.2) vanish, Ike tran^om 

will take i^s simplest form and will ue reduced 

multiplication of initially choeen solutions by certain 

constants, ihese are sometimes hnown as Ploquet solutions. 


(A-2.3.1) 
(A— 2.3 .2) 


y*(t+'l) = P x 

y*(t +'1) = p 2 y* 2 (t) 

prom the theory of linear transformations, the trans- 
formation of the type U-2.2) can be reduced to the diagonal 
form (A-2.3), where the constants P lj2 are known from 

characteristic equation 


a il ' P 


a. 


•21 


12 


a 22 P 


0 


whi ch can 


he written in the following form, 


2 - 2Ap + B = 0, 


(A-2.4)- 
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where, 

A = I ^ a ll + *22*’ 

, B = a ll a 22 " a i2 a 2l ’ 

IT he characteristic equation (A-2.4) defines the 
character of the solution of the equation U-2.1) . Hence it 
is desired to get the form of equation (A-2.4) from the 

independent linear solutions, y 1 ( t ) arL ^ * Let US 

assume the most general initial conditions xOr 

y ( t ) as follows. 


y x (o) = 

y-i(°) = 


y 2 (0) = a 2 


J 2 (o) 


( A-2.5 ’ 


where (■) denotes differentiation with respect to t. i'hes.- 
initial conditions ere required to be satisfied by the 
solutions (A-2.2). Hence one gets using (A-2.2) and (A-2-5) 


a ll “l + a 12 “2 = 

a ll h + a 12 ? 2 = h (I) ’ 

a 21 a x + s 22 a 2 = y 2 (I) ’ 

a 21 h + a 22 p 2 = y 2^ ‘ 


(A-2.b .1) 
( A— 2 . b . 2 ) 
(A-2.6.3) 
(A-2.6.4) 


Equat 


ions (A-2.6) determine the constants a ±; 


follows . 



154 


a 2 ^1^ ~ ^2 y l^ - 

a ll a 2^1 — a l^ 2 

i 3 -^y-]_( i) ~ a 3 _y]_( x ) 

a 12 = " a 2 p 1 

a p y ? (g) - p 2 y 2 (i) 
a 21 = a 2 p x - a x p 2 

Pl y ? (l) - a^d*) 

and a 22 - 

_ ~ ^-tsnts a , the constant B in (A-2.4) 

Using the values of constants a j_-j 5 

is given by 


B = a n a 2 2 " a !2 a 21 


y,(!) y 9 (i) - y 2 (?) 


•'2'~ / " 2 
;«iP 2 - a 2^r 


( A-2.7 ) 


in the following it is proved that B is always equal 
(t) are solutions of equation (ii-s.l) 


to unity. Since y ly2 
hence, 


*1 + 


q ^ [i - 2s f(t) ]y x - 


y 9 + ^ 


o 2 [i - 2e f(t)]y 2 


(A-2.8.1) 

(A-2.8.2) 


Multiplying (A-2.8.1) oy y 2 

from the other, one gets 


(t) and (A-2.8.2) by y-^t), 


and subtracting one 
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y x (t) y 2 (t) - y 2 (t) h (t> = 
On integrating once, one obtains, 
yi (t) y 2 (t) - y 2 (t) • 


C A— 2-9) 


where C is the integration 


constant 


Using init 


ial conditions (A-2.5j.( at ^ 


tj; 


on 6 obtains 


C = a-$ 2 ~ a 2^1 ’ 

and at t= r i , equation (A-2.9) g iveS 

yi (D j- 2 d') - y 2 (x) f 2 (1-) = e 

Hence from (A-2.7) > 


= a x p 2 - oc 2 P x . 


B = 1* 

characteristic equation U-2.4) becomes 


:‘hus, the 

p ^ - 2Ap + 1 = ^ 


(A-2.10) 


and 


its roots are relate^ 


P 1 P 2 


(A-2.ll) 

the particular 


It has been alreaay shown tnas enong 

„f (t-2 1) two linearly independent solutions 
solutions o± vA-^.r; 


y * (t) exist wtn 
-L * 2 


which satisfy equations 


A-2.3) • 



yj(t«) = 4 (t) * k - 1,2 


These solutions, which acquire a. constant multiplier oy 
the addition of the period 1- to t, can be represented in 

the form 


(x) «x„(t> . *=1,2. 


(A. 2.12) 


where Xl _ 2 (t) are oertsin peri0dlC tmci±ona ° f Pen ° d X '' 

It follows from equation U-2.12) that the behavior 

ox the solutions as depends on the value of charac- 

tV _ vroiue of its moduli). 

teristic roots (more precisely? 

In fact, taking into account tnat 

in p = In | p 1 + i s.rgp ? 

f a p 12 ) in the following form 
ie can write equation (A-2.12) 

( t /T) In | Pk ! k = i f 2 (1-2.13) 


on 


y k (t) 


fl k (t) e- 


where „ 


, , , + i , (it A 1 ) atgp 

f lk (t) = X k (t) = 


and is bounded. 

, . 0 1H \ on , can establish the stability- 

looking to ( ii.-2.lb) on-- call 

Tf the characteristic numDer p k >] -* then tn,, 
criteria. H biiC * 

corresponding solution U-2.13) Deeses unboundedly with 
time. If p k < 1, then a bounded decaying solution wxtn ti„e 
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is obtained. If P k = 1, then a periodic and bounded sola- 
tion with time is obtained. 

B is always equal to unity as noted earlier. Hence 
on uhe basis of the characteristic value p k one can establish 
the limits on the values of A for stability using equation 

( A-2.4) • 

If j k\ >1, then the characteristic roots will oe 
real, and one of them, will be greater than unity. Hence 
the general solution of equation (A-2.1) will increase 
unboundedly with time. 

If | A} <1, the characteristic equation has complex 
conjugate roots and since their products must be equal to 
B=l, their moduli will be equal to uni oy (nq. A 2-11) 
the case of complex characteristic roots corresponds to the 
region of bounded solutions. On the boundaries separating 
the regions of bounded solutions, condition Ui = 1 must 
be satisfied [Eef . [13] ] .- 

Por | = I, multiple roots occur, moreover, it 

. . / , n - j - \ thot cnch roots can be 

follows from equation (A-2.1^J tha. snn xu a. 


either = P 2 ~ ~ orp l 


o = -1. In the first case 
m 2 


as seer from U-2.5), solution -ill b« periodio with 

a period i = |fi whereas in she second case the period 

will be 2i • 
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therefore, the regions of unbounded solutions are 
separated from the regions of bounded solutions by the 
periodic solutions with period i end 2x- I" other words, 

„o solutions of identical period bound the region of 
instability ( |A|> 1), and two solutions of different 
periods bound the region of stability (| A|<1)- 

Bolotin [13] has shown similarly for the coupled 
pype of homogeneous Bathieu-Hill equations, that the solu- 
tions of an identical period bound the region of instability 

„ o . r fcTt(Stl t periods bound the region of 
and tw solutions of diiferent peaoc 

stability . 
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APPENDIX B 


30LIP1T0N OP MA’l'HISU'S EQuMilON 


Following Mat hi eu ’ s equation is considered to find 
the solution by small perturbation met nod. 


x + Q 2 (1 - e cosmt) x 3 A. 


( B.l) 


in which A is constant and e is a small parameter. Initial 
conditions for (B.l) are 


x(0) = x(0) = 0 


(B.2) 


One assumes the solution of (B.l) in a power series of e 
as given below. 


x = x o +ex-^ + e x^ + 


-r £ k x. 


(B.3) 


On substituting (B.3) into (B.l) and (B.2) and comparing 
the powers of e k on both siaes of tr.e equations, one obtains 


x^ + a x = 

o ° 


with, x(0) = x(0) = 0, 


(B.4.1) 
(B.4. 2 ) 
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and 


with, 


x k + ft 2 x k = ft 2 X k _ x cos u t, 


X, 


Xo) = xjo) = 0 , k = 1 , 2 , 3 , . .. 


(B.5.1) 


(B.5.2) 


k k v ~' “ “k v 

The solution of (B.4.1) with initial conditions (d.4-2) is 

x (t) = ~ (1 - cos at) (B * 6) 

O ft 2 

Similarly the solution of x L using (B.6) one obtains from 
(B.5.1) and (B.5-2) , 


cos u t . 1 , cos(^.+m)t, _ c os(,o. -o) t ) 

o o » o V- i > 1 ii\ { PD — m) 


z , N . A r cub op ^ ± f v • 1 , ..y- - —rrr — r 

(t) = A [ o 2\ + 2 ^ ( 2fi+to)W o l Hi -w) 


(f2 -a) ) 


d£r cos nt 


(n 2 -m 2 )(4o 2 -w 2 ) 


] 


(S.7) 


Proceeding the same way as for x L 


one obtains the solution 


of Xg as 


,2 


x P (t) = *§- l ~ - 2 2w" 2 ~ > ^ 

2 2 (Q -a> )(4S1 -to ) 


: - OsCn -to)t _ cos (Q +CJ ) t ',\ 
v "7 O.-.L.. . 'l ( 2Q + ov ^ 




1 / cos(^ +2m)t , cos(Q -2 apt \ + — 00 p 2 -^ d 

~ “2 VC2n+toililW (20-coKS-eol' , Q 2 _ 0 , 2 )( o 2 -4a> 2 ) 


(Q 2 - 9tofl 


+ 


0) 


2 ( 4fi 2 -a) 2 ) (fi 2 -4<o 2 ) 


- COEfit ] 


(B.8) 


One may confine up to the terms containing e* 
of x and thus approximating x as , 


in the series 
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mmm 2 -1 

x + ex 1 + e x 2 


(B.9) 


. « y x in (B.9) ons obtains 

On substituting for x 0 * x 2 1 v 

i > / -3f£cosJ_t___ cos jot 

x = A [ij U - cosnt) + < (a 2. u 2 )( 4 c 2^2 ) (a 2 . u 2 ) 


1 ,cos( a -w)t cost a- M)t s \ x r 2 — ( — 5 — 57? ? 27 

+ 2 M20+1) u ( 22^7) V 2 -U )(4a -a ) 


cos (fl- Kd)t _ cos( 0- ^ V ) + 

( ,,U oO jji! 1 W ( 20-OJ) 




003(0-2011^)+ cos 2 i ^ 

(2Q-03TTO-03) (q 2 - w 2 )(2 -4 o> ) 


( n 2 -9a) 2 ) cosOt 

u 2 (40 2 -u 2 )(0 2 -4w 2 ) 


(B.10) 



